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O círculo de precisão rotring 


COM COMPASSOS DE PRECISÃO ROTRING 
Caso pretenda traçar círculos de precisão, então necessita dum compasso 
de precisão. Dum compasso extraordinário. Dum compasso facilmente 
regulável e que não oscile. Com um comando de simetria exacto. Para 

um trabalho sem esticões. Para círculos mínimos até 0,7 mm de 
diâmetro — neste caso recomendamos o compasso de círculos 

mínimos RP 06. E para círculos grandes até 530 mm de a 
diâmetro. Agora, em poucos segundos pode passar-se da 
mina para a tinta-da-china ou vice-versa. 
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desengatar o parafuso micro- 
métrico, fixar a perna do 
compasso na posição 
desejada — dar um 

ajuste ao parafuso micro- 
métrico — e o compasso 

RP 01 está bem ajus- 

tado. A superfície deste 


extraordinário compasso é 
é cromada e muito : 
brilhante. 


Art. No. 532 101 
(com alongador, 
porta-lápis 

e pontas secas) 


“2. COMPASSO TÉCNICO RT 04 
ROTRING 

Um compasso com peça intermutável 
em forma de U, em que se coloca 

a ponteira de tinta-da-china, podendo 
esta, no entanto, mover-se. Dois 
segmentos dentados, encaixando 
exactamente um no outro, garantem 


Na perna da ponta seca encontra-se um 
alongador telescópio. À superfície do 
compasso é cromada e baça. 

Art. No. 531 104 
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eixo central. Com este compasso pode 
traçar com toda a precisão e limpeza 
rotring-werke Riepe KG, D-2 Hamburg 54, P.0. Box 541060 círculos mínimos até 0,7 mm de diá- 
metro devido à não-oscilacão da forte 
mola de lâmina. 
A superfície do compasso é cromada 
e de elevado brilho. 
Art. No. 532 106 
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Do conhecimento e da acção 


RESUMO 


O modo como o homem adquire conhecimento do 
universo, como o conceptualisa e actua. 


1 — INTRODUÇÃO 
1.1 — O UNIVERSO DOS SERES E O DOS ATRIBUTOS 


Um ser (S) (conceito primário não definido) possui 
vários atributos (A) (conceito primário igualmente não 
definido). 

Os seres pertencem ao Universo dos seres Us 

Os atributos pertencem a outro Universo, o dos atri- 
butos U A: 

O homem conhece os seres pelos seus artibutos. 

Os seres, elementos de Us projectam imagens no 
universo dos atributos U, por meio de Relações que 
simbolizaremos por p com índices apropriados. 

Os seres interaccionam uns com os outros e dessas 
interacções resulta que as respectivas «imagens atribu- 
tivas» são alteradas. 

Postula-se que dois seres não interaccionam se as 
respectivas imagens no universo dos atributos não se 
modificam como resultado dessa interacção tentada. 

Figurativamente podem apresentar-se estes concei- 
tos da forma seguinte: 
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ANTÓNIO GOUVEA PORTELA 


SUMARY 


How man aquires knowladge about the universe, 
conceptualises it and acts. 


Na Fig. 1 está representado um ser S e um conjunto 
de Atributos A , que S projecta por meio dum conjunto 


de relações c; 


Fig. 2 


Na Fig. 2 estão representados dois seres S, esS, 
e as respectivas constelações de atributos A ,; " As; 
obtidos por meio das duas constelações de relações 
o dh 

Se entre S;, e S; se se estabelecer uma interacção 
então, pelo menos, uma imagem atributiva A:;, e outra 
As; se modifica o que demonstra que uma interacção 
tevo lugar, como se representa na Fig. 3. 


(O Re + 21 | Y <) y 
Ao RS 
MT ; Go) Tem Lai, - Cáreis da 
a do Zi a 
F, | / | À | | 
4 k Aa 
f, / o 
Fig. 3 


1.2 — INSTRUMENTOS (SERES INSTRUMENTAIS S , ) 


Um instrumento (no sentido lato) não é mais que 
um ser que tem uma imagem atributiva «típica». 

Um termómetro é construído para medir o artributo 
da temperatura. 

Um instrumento é optimisado no sentido de dar 


relevancia à imagem que projecta num certo atributo Aa . 

Procura-se ainda que, quando usado, isto é, estabe- 
lecida a interacção do ser-instrumento com o ser-mensu- 
rado, cujo atributo se deseja apreciar, se verifique que 
os atributos do ser-mensurado não sofram modificações 
relevantes de modo a que se possa declarar que este não 
foi submetido a uma alteração por força da interacção, 
com o ser-instrumental. 

Assim, o homem foi conduzido a fabricar uma vas- 
ta classe de instrumentos cada um especializado na 
apreciação de certo atributo. 

Figurativamente a imagem-atributiva de um ser não 
é, na ciência moderna, conhecida senão por intermédio 
das imagens típicas dos seres instrumentais que são 
postos em interacção com ele. 
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Fig. 4 


Na Fig. 4 estão indicados os passos normalmente 
dados na operação da observação-atributiva de um ser S. 


NOTA 1 -— O homem pode desempenhar o papel de ser- 
-instrumento, quando propõe um juizo de valor. Deste 
modo o homem desempenha duas funções, a de Ser- 
“instrumento e a de apreciador de atributos. 


NOTA 2 —- Este conceito é muito vasto e se é simples 
reconhecer que um termómetro é um ser-instrumento, é 
mais difícil reconhecer que uma simples contagem dos 
elementos dum conjunto ou a verificação de que a luz 
está acessa, são para este efeito, também instrumentos. 


NOTA 3 — Sem querer introduzir qualquer conceito de 
«informação», pode dizer-se duma forma figurada que na 
interacção entre S, es iK Se podia desprezar os efei- 


tos | «o € suprimir no grafo os arcos desse tipo. 


pa 


Fig. 5 


Num instrumento procura-se que o grafo da Fig. 5 a 
seja substituído pelo grafo da Fig. 5 b. 


1.3 — SERES INFLUENTES: (OU POTENTES) 


Em oposição aos seres-instrumentais que se carac- 
terizam por evitar, na medida do possível, alterar O ser- 
-mensurado quando aqueles são postos em interacção 
com estes, outros há que se destinam justamente a pro- 
vocar modificações nos seres-mensurados. 

Com efeito espera-se que uma válvula, um contacto- 
-eléctrico, etc., interfiram no ser em estudo de modo a 
alterar a sua imagem atributiva. 
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Fig. 6 


Um Ser-Potente S p interacciona com um Ser dado 
Ss as Então deve verificar-se que: 


- A, imagem-atributiva do Ser-instrumental S ;, 


em interacção com Sp 

FA, imagem-atributiva do Ser-instrumental St em 
interacção com 8, 
se modifiquem por força de interacção, entre S pº So 

Na Fig. 6 está representado este conceito e pro- 
curou-se ainda distingir os S, destinados a obter ima- 
gens-atributivas de Sp e S, e cuja intervenção pouco 
devo alterar a interacção em estudo de Sp com Sa: 

Note-se finalmente que no grafo que descreve a in- 
teracção entre S, e S, tanto RA como io são em 
geral importantes. 

Contudo pode acontecer que Ver possa ser despre- 
zável, (ao contrário do que sucedia com |, ). 

Assim, teremos as seguintes situações típicas: 


—— E, 


Loo Lor 


Fig 6a 
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Na Fig. 6a, está representado a situação real, isto 
é, a interacção envolve sempre um duplo arco. 


— 


4oT 


LPo 
(se YE) 


Tor 


Fig. 6b 


Na Fig. 6 b, para instrumentos é possível desprezar 
na generalidade o arco | ,,. 


Na Fig. 6 c, para uma relativamente vasta gama de 
Seres-potentes pode desprezar-se | , O grafo obtido 
constitui uma cadeia linear. 


op 


| 


DE 


1.4 — RELAÇÃO ENTRE IMAGEM-ATRIBUTIVAS 


Fig. 6c 


Há três grandes categorias de relações entre ima- 
gens-atributivas. 


a) Categoria f. 


Seja dado um ser S e suas respectivas imagens- 
atributivas A, , ie Q 


Se for possível definir uma relação f tal que 


f 
di Se M 


onde: jeQ,cN eo leM, 
onde 2, é o complemento de 2, em relação a $., pode 
então declarar-se queas imagens A, estão entre si rela- 
cionadas, (ou correlacionadas). 

Esta categoria de relações só demonstra a existência 
de redundância na informação colhida a respeito de S e 
têm interesse sobretudo para controlar (check) a informa- 


ção colhida. 
b) Categoria F. 


Relações que têm por domínio Imagens-atributivas 
de um ser e por contradomínio outros atributos (em geral 
não mensuráveis). 

Exemplos: energia, informação, entropia, liquidez, 
prezer, tensão nervosa, etc. 

Este tipo de relações têm a maior importância 
porque permitem descrever os seres reais à luz de atri- 
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butos conceptuais, isto é, do foro dos conceitos e não 
dos sentidos. 


c) Categoria q 


Relações que se estabelecem entre imagens de seres 
interaccionando entre si, nomeadamente entre seres po- 
tentes e o ser em observação. 

Tem interesse figurar estes três tipos de relações. 


Fig. 7a 


Categoria f 


Seja aoj* A; com je £2%, o valor do atributo de S 
em A,. 


o 


Uma relação do tipo f, terá a forma 


f 
Ao do % Àoj 


com: jeMcNoeleg, 
Como por exemplo a, = 2 ao 
jeS, 

então ã,j não traz nova informação e o conhecimento a oi 


je 92,era suficiente, a, é informação redundante, que 
pode apenas servir de controle. 
X A j 


O Domínio de f é o produto cartecano| je 
| 


A 
e Contradomínio é | 
es2, 
veja-se Fig. 7a 
Categoria F 
Dom F = X A j onde 2, e QM 
je, 


e cont F =X AB 
Be ç 
Porque as imagens-atributivas Af com f «(% são as 
mais variadas pode até acontecer que se desprezem 


completamente os A e se substituam por AR para 
tanto basta que F seja invertível. 


Us 
A 


e aa RIMAS 
ip pç 


1 mia 4 Lu 
S Lex da 
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Fig 7b 


Na Fig. 7b, está representada a forma como se 
criam estas imagens-atributivas secundárias a partir das 
imagens-primárias. 


Fig. 7 c 


Categoria «q 
— 55 interacciona com S ,, 
— projecta imagens do tipo A, com ves 


e S Projecta imagens do tipo A, com Ke6 


tem por domínio X Ar 


re meu 
Y 
e contradomínio X A |. 
k 2 6,co 
a 
tem por domínio X A 
k e 63C À 
Pa 


e contradomínio X Ar 
Fr ées,c su 
Estas relações q, no universo Us das Imagens- 
atributivas, procuram descrever, tanto quanto possível, 


O que se passa entre S pº S no mundo do Real. Veja-se 
Fig. 7 c. 

Convém ainda chamara a tenção para os seguintes 
aspectos: 


— So designarmos por S, à reunião interactuante 
de SUS p S y . então as relações do tipo 
q confundem-se com as do tipo f. 
Sobretudo, em física e automação é costume re- 
correr alternadamente a estes dois modos de en- 
carar o par (S, S p)- 


— So designarmos por AM à reunião interactuante 
dos seres IS, ) US, USUIS, 1 -S, 
e, por outro lado reunirmos |A;) U IA ] A 
então A,éa imagem de S por | Ep ) com 
| v 6 u !; Veja-se Fig. 7d. 


Se na resolução de um determinado prob- 
lema a imagem pw fz, | for adequada então 
'A w)] deverá conter toda a informação necessá- 
ria para resolver o problema proposto. 


e 


Fig. 7d 


— O estudo dum fenómeno interactivo poderá fazer- 
-s8 como se indica na Fig. 7 e, o que já é uma 
representação mais simples. 


um 
a, 


Fig. 7e 
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— Por último tem interesse referir que o formalista 
desinteressa-se do universo Uc, alguém (físico, 
engenheiro, economista, médico, astrónomo, etc.) 
so encarregou da tarefa de formar as imagens 
(A w). O seu problema nasce daqui por diante. 
No capítulo seguinte 1.4) trata-se sumaria- 

mente deste problema. 


1.4 — FUNÇÕES DO FORMALISTA 


No contexto desta exposição o maior problema do 
formalista é tornar nítidas as imagens (A) o que através 
do tempo foi obtido usando várias linhas de progresso: 


A — Definição sintática de cada um desses espaços 
onde a imagem A, é descrita, — aqui temos 
toda a álgebra e topologia. 


B — Oferecer ao físico (muitas vezes por iniciativa 
do próprio físico) uma grande cópia de espa- 
ços: grupos, anéis, módulos, espaços funcio- 
nais, etc., com métricas variadas (e sem elas), 
conceitos de distância, distribuição, etc. 


C — Estabelecer regras de inferência de aplicação 
sintática: (operações, operadores) que dum 
modo geral chamaremos de Relações, que ha- 
bilitem a racionalizar o juizo sobre as imagens. 

Qual o principal perigo? o confundir-se a Nuvem 

com Juno, isto é, a Imagem do ser, com o próprio ser. 

Leis da física, modelos econométricos, caixas pretas 

da álgebra de Blocos, são Relações, entre imagens de 
interacções reais de entes reais, em geral purificadas, 
idealisadas e sempre com um conteúdo de informação 
inferior ao do Real. 

O Estudo deste problema é o objecto do capítulo 

seguinte. 


2 — DO CONHECIMENTO E DO CONTROLO 


O Homem presistentemente persegue o objectivo de 
Controlar os acontecimentos. 

Para a consecução deste objectivo é necessário 
conhecer as variáveis potentes ou influentes sobre as 
quais pode actuar. 

Para atingir esta segunda finalidade, que aliás nem 
sempre consegue, é mister resolver duas classes de pro- 
blemas: 


a) Primeiro — Adquirir «conhecimento» a respeito 
do ser ou fenómeno porque se interessa. Escolher o es- 
paço atributivo adequado. 


b) Segundo — Adquirir «conhecimento» sobre as 
variáveis influentes no ser ou fenómeno e sobre as variá- 
veis que pode manusear (comandar) que directa ou in- 
directamente interferem nas variáveis influentes referidas 
em 2º, lugar. 


a) Primeiro Problema: 


Foi tratado no capítulo anterior, e já vimos como 
é possível criar um espaço imagem-atributiva para ca- 
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racterizar um ser (um fenómeno ou interacção) conve- 
nientemente. 


b) Quanto ao Segundo Problema vamos passar a 
expôr: 


Seja Sy um ser (ou conjunto de seres eventual- 
mente interaccionantes) e seja (A | K «6 a sua ima- 
gem atributiva obtida por meio de relações 2 + com k e 
9 e seja tA,, |). v e 2 a imagem dos seres influentes ou 
potentes por ( c,),r e 2; Sr inclui os seres instrumen- 
tais necessários, bem como os seres Potentes e respec- 
tivos seres instrumentais. Na Fig. 8, representa-se o 
mesmo que em Fig. 7 c, mas engloba-se em S y tudo 
quanto está representado e analisado em U da Fig, 7 c. 


Fig 8 


Escolhamos os espaços formais adequados para re- 
presentar as imagens-atributivas: 


oa K ec 6 
X + nd res? 


e formaremos o respectivo produto carteziano: 


Keg re 


o finalmente, seja q um conjunto pertencente ao Produto 


Carteziano seguinte: 
3 X ú I T X ” 
7“ reM keg 


y relaciona os dois conjuntos multiplicados carteziamente 
Emgeral sucede que não é possível criar, duma só 


” X 
vez reQ , ) + porque só a pouco e pouco é possível 


adquirir o «conhecimento» de quais os atributos que efi- 
cazmente são influentes; estamos no domínio da ciência 
experimental. 

Também não será dum só jacto que se saberá es- 
colher convenientemente a estrutura algébrica e topoló- 
gica tanto para os espaços X, como para o espaço 
Te X+ 
rel? 


Problema igual, mas em geral mais simples se veri- 
fica para  X 
Ke 
Neste ponto o que se deseja realçar é o modo como 
progressivamente esta operação pode ser efectuada. Com 
efeito, supunhamos que= X, já foi identificado e es- 
Keg 


ni 


colhido. 

O Problema está em que foi escolhido um conjunto 
1), de coordenadas e ao tentar-se descrever a função q: 
obteve-se uma relação com a forma indicada em Fig. 9 a 
(Figura simbólica). 


TT xy 
re, 
Fig 9a 


O Controle do fenómeno (ou Sêr) é fraco e prova- 
velmente não satisfaz às necessidades, com o efeito o 
sistema pode sofrer grandes alterações para o mesmo 
conjunto de Comandos. 

É essencial procurara mais variáveis influentes e 
Q,=89 é o novo conjunto escolhido. 


xr 
FE 


Fig 9b 


Na Fig. 9 b está representada a nova relação q, é 
verifica-se que a região de incerteza diminui substancial- 
mente. 

Um terceiro conjunto 2; > 4 o $S! foi adoptado e a 
relação obtida q; já será suficientemente boa para não 
prosseguir nesta operação: aumento do conjunto 2, de 
variáveis influentes (Fig. 9 c). 


Tx 


KEO 


Fig. 9 c 


Tem interesse referir os motivos que levam a não 
prosseguir a operação: 


1 — Dificuldade de encontrar novas variáveis in- 
fluentes ou potentes. 
Com efeito aquelas cuja influência é mais 
evidente e relevante foram as experimentadas 
em primeiro lugar. 


2 —. Se o número de variáveis influentes for muito 
grande o seu respectivo controle torna-se pro- 
gressivamente mais difícil. 


3 — As vantagens marginais que advêem dum con- 
trole mais correcto da imagem de Spemz X, 
rel) 
são cada vez menores, 


c) O Terceiro Problema: 


O facto de se explicar a imagem de Sy por meio de 


uma outra imagem em= X,, não significa que se possa 


Fe “2 
controlar a imagem de Sr. 
Exempos óbvios são situações tais como: o movi- 
mento de Rotação da Terra, a atração lunar, etc. 
Assim muitos problemas deixam de ter solução so- 
mente porque o homem pode apenas explicá-los mas não 
controla-los. 


Porém há casos menos fortes do que os acima men- 
cionados, e que se situam no interior da sua fronteira de 
acção e que podem formular-se do seguinte modo. 


Conhecem-se as variáveis influentes = X pr mas não 


Keg 
se sabe controlar todas ou algumas variáveis, então há 
que procurar novas variáveis que possam substituir as 
primeiras mas que sejam já susceptíveis de ser contro- 
ladas. 

Na Fig. 10 figura-se entre real SG que interaccio- 
nam com Sy igualmente real, de tal modo que as ima- 
gens em (A, ) e (Ag), em conjunto, não só explicam a 
imagem em (A) como ainda, algumas delas corres- 
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Fig. 10 


pondem a variáveis potentes que eficazmente controlam a 
imagem em (A,) de Sp . 
Se designarmos por Za os elementos pertencentes 
que: a Xy como a X, onde Za é potente e daí que pode 
ge É req 
ser controlado, teremos que haverá um conjunto F tal 
que q Xa é um espaço suficiente de variáveis potentes 


ae E 
(Z.: X, ) que permitem controlar a imagem de Sr em 
TX. 

K 
Keg 


Note-se que Card E > Card S), com efeito, se tal não 
sucedesse, a qualidade do controle seria prejudicada, 

Exige-se, em geral, ainda que exista uma relação 
entre a imagem em 7 X, de Sy e a imagem de (Sg USy) 

re!) 

em 7 X q, do tipo bijectivo. 

g ef 

Esta segunda condição garante que entre as variá- 
veis potentes e as variáveis explicativas existe uma rela- 
ção bijectiva, 


3 — RELAÇÕES 


Nos capítulos anteriores ficou patente o papel iIm- 
portante que as relações desempenham na manipulação 
da informação. 

Sob o ponto de vista informação e controle convém 
classificar as relações em duas grandes classes: 


as que destroiem informação 
as que preservam informação 


Além desta classificação haverá todas aquelas, que 
são correntes nas teorias formais como a dos conjuntos: 
aplicações, funções, medidas, etc., etc., que supomos 
conhecidas e portanto só nos vamos referir à classifica- 
ção acima por ser esta que mais importância tem para 
o tema em exame, 


3.1 — RELAÇÕES QUE DESTROIEM INFORMAÇÃO 


Comecemos por dar um conjunto de exemplos para 
ilustrar este conceito. 


— Os seres-instrumentais, uma balança, um termó- 
metro, etc., fornecem por interacção com um ser 
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-mensurado dado, uma informação típica; um 
peso, uma temperatura,etc. 

O instrumento e o respectivo modus-ope- 
randi, consideram-se tanto mais afinados e efi- 
cazes quanto menos perturbada for, a informa- 
ção, pelos restantes e infinitos atributos do ser- 
-mensurado. 


Tudo se passa como se toda a restante in- 
formação a respeito dos atributos do ser-mensu- 
rado fosse destruída (ou perdida) para só ficar 
aquela para a qual o instrumento foi concebido, 


— Um filtro monocromático, só deixa passar a ener- 
gia associada à frequência para o qual foi cons- 
truído. 


— Reduzir a informação a respeito de um conjunto 
de seres, apenas ao conhecimento do seu nú- 
mero é também destruir informação. 

Também há inúmeros exemplos no domínio do uni- 

verso dos atributos (U A). 


— Escolher um vector extremando, num programa 
linear, elegendo uma funcional para critério único 
de escolha. Toda a informação contida nos ele- 
mentos dum vector é reduzida a um real. 


— O Produto ou a soma de vários números, reduz 
a informação, com efeito, há mais informação 
conhecendo todas as parcelas ou termos do que 
apenas o seu produto ou soma. 


— Operações do Cálculo de afinores, tais como, 
transvecção, contracção, produto misto ou alter- 
nado, destroiem informação. 


— À projecção, o corte, a sombra de uma forma 
geométrica têm menos informação do que a da 
referida forma. 

Dados estes exemplos tem interesse explorar um 

pouco mais dois casos: 


Funcionais 
Distribuições 


3.1.1 — FUNCIONAIS (OU FUNÇÕES GENERALISADAS) 


Seja dado um espaço produto cartaziano x, Xx: e 
seja R uma relação, Rox, Xx: e representemos por 


C, Oo corpo dos Reais. 


Pode conceber-se uma Relação F que tenha por do- 
mínio o conjunto assim descrito: 


dom =! (x;, Xj ) & R) 
| 


Isto é, todos os pares (x; | X; ) que pertencem a 
R e têm por contradomínio os reais. 


Então F reduziu toda a informação contida na Re- 
lação R a um real a. 


TAM” 


Esta operação pode ser usada de dois modos: 


a) ou Procurar todos os R, que podem constituir 
dom F e formar com eles um conjunto: 
UR, :R; é dom F) 

Esta operação é típica para as funcionais de pro- 
gramação linear ou outra. 


b) ou Criar um conjunto de relações F. cujos ele- 
mentos têm entre si certas familiaridades e a partir dum 
mesmo R, obter os a ii correspondentes. 

Ora como cada «jj está associado ao índice j que 
caracteriza a relação F, poderá dizer-se que desta ope- 
ração resultou um par (aj. j). 
onde | cão ú : 


jeQ 


(£2 é frequentemente o corpo dos Reais C, ). 


e então R, [CN (ao DM) 


isto é, R; eo conjunto de relações | f, | deram 
origem a uma relação OQ, cC, X “2 


tal que 


agé a imagem de R, por F. 
Qi = (args j): R s 
jé o índice do F ; 


Pode então constituir-se um conjunto de relações 
tF,) suficientemente «denso» de modo que toda a in- 
formação contida na relação Ri seja transferida para a 
relação Qi, tornando possível assim inverter-se a opera- 
ção de modo a regressar a Ri a partir de Qi. 


Neste caso, teriam sido construídas Relações 
(F;) e tF, ) = que não perdem informação. 
c) ou Finalmente, o par de conjuntos de relações 


F,) e (F, podem servir para tratar conjuntos de 
relações (Ri), (aos quais fossem aplicáveis), e produzir 
as Relações do tipo (Qi) e inversamente, 

Como exemplos temos as transformadas de Fourier, 
de Laplace, os momentos, etc. 


3.1.2 — DISTRIBUIÇÕES (PROBABILIDADES E ESTA- 
TISTICA) 


a) — Regressaremos ao 2.º capítulo e ao problema 
tratado na alínea b) e em particular às figuras 9-a, 9-b e 
9-c. 

Nem sempre é possível reduzir a região da incer- 
teza pelo processo de acrescentar mais variáveis X, 
explicativas (ou potentes). 

Pode até acontecer que a relação seja de tal modo 
fraca que se confunda com o produto carteziano do seu 
domínio pelo seu contradomínio, então é possível às 
vezes (quando o contradomínio de q é probabilizável) 
fazer corresponder a cada elemento (...x, ...) e mx 


Fr 
keg 
não um elemento (...x,...) eq x, mas uma relação 


r 
kea 


r 


Pre assim construída. 


Domínio P, em x 
K 
keg 
Contradomínio P, um real p, e R+ 


Se esta relação P. fôr do tipo medida normada com 
norma igual 1, pois temos construído o tipo das distri- 
buições usadas em Probabilidades. 


Fig. 11 


Na fig. 11 está representado o caso de uma só va- 
riável explicativa (ou potente) x, e uma só variável a 
explicar Xe Assim temos: 


— X, (variável explicativa) que na Fig. 11 se repre- 
sentaram os valores x,| e x ambos perten- 
centes a X + 


r2 


A (variável a esplicar) 
— p uma variável real positiva, e p e [0,1]. 


—. Pei [3 Pros duas relações (X (o p ) que estão em 

correspondência respectivamente com X 8 Xo 

Estas relações P., e P,, têm um conteúdo de infor- 
mação importante, pois pode conhecer-se qual o valor 
X 4 mais provável para um certo X + dado ou a média 
ou a respectiva variância, enviezamento, etc., etc. 

Para o efeito basta aplicar certas funcionais (os 
momentos) para fazer a correspondência entre cada uma 
das distribuições P,; com um certo x : Reais. 

O Processo completo pode, como exemplo, aplicar- 
-se ao valor médio. Veja-se a Fig. 12. 


Xu «R 


Fig. 12 


Em muitas aplicações, quando a variância é aceitá- 
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vel, pode confundir-se a aplicação sucessiva das duas 
relações (y o À) por um só p: 


a e NE SEO 


No exemplo apresentado, x, admitiu-se dado o de- 
terminado e x, aleatório. 


onde 


b) Tem interesse comparar o conteúdo de informa- 
ção entre uma relação com uma incerteza irreductivel e 
a informação adicional que pode ser carreada se esse 
domínio de incerteza poder ser probabilisado. 

Novamente usaremos na ilustração um só x. e um 
SÓ Xp 

Tratado X, como variável determinada, a relação «w 
toma o cspecto da Fig. 13 a. 


Xv 


Xy 

Fig 13 a 
É claro que x, explica (ou comanda) mal x, . 
Tendo sido possível probalisar o espaço assim 


definido: 


Xp Xké 2.º membro de q quando x, é 1.º membro) 


“ 


Fig 13 b 


As Distribuições respectivas Fig. 13 b contêm infor- 
mações suficientes para declarar, por exemplo: que com 
uma probabilidade de 90 % a Relação «q entre Xx.º x, é 
a da Fig. 13c. 
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X 


D = 40% 


Xw 


aô Fig. 13c 


Kw 


40% 
19% 
I00,/ 
Xe 
Fig. 13 d 


Na Fig. 13 d estão representados os contornos que 
delimitam as regiões dos espaços x, X X onde q está 
em correspondência com as probabilidades referidas na 
dita Fig. 13 d, isto é 90 %, 99 %, 100 %. 


c) Finalmente tem interesse referir à Corelação en- 
tre dois XX 48 X 


Xe Xo são duas variáveis «insuficientemente» 


X 
res) 
Pergunta-so x,: e Xp2 São ou não independentes 
entre si ou existe um certo grau de dependência? 
So X,1 6 Xp: fossem determinadas (ou tratadas 
como tal), então a relação: tomaria a forma dada na 
Fig. 14 a e o contorno indicado poderia no caso extremo 


explicadas por 


2 
confundir-se com 7 Xi 
jel 


Xw 


Xuz 
Fig 148 


t 


Pode porém ser possíve: estabelecer-se em 7 X; 


lel 
uma medida normada (Probabilizar o espaço referido) e 


então obter-se-ia, por exemplo, a Fig. 14 b. Assim alguns 
pares (Xu; e Xi ) estariam associados a valores de p 
maiores do que outros, e novamente seria possível 
extrair maior informação da Fig. 14 b, do que da Fig. 


14 a. 
Xu 


Fig. 14b Ku 
RE 


Se designarmos por:' a relação: X a X X,o e* p 
então o" permite extrair informação adicional. 


+, 


Resumindo o Capítulo 3.1: 


a) As relações que destróiem informação são essen- 
ciais para o refinamento do nosso conhecimento a res- 
peito do real. 

São as que permitem extrair formas atributivas 
puras. 

O universo dos atributos foi construído a partir do 
universo dos seres por meio de relações (em geral em 
hard-ware) que têm esta propriedade filtrante, desenvol- 
vida em alto grau. 

Já no domínio dos Atributos, também são as rela- 
ções deste tipo que permitem escolher e ordenar. 


b) Mostrou-se como as funcionais têm esta proprie- 
dade em alto grau e ainda como é possível que conjuntos 
de funcionais pudessem em certas circunstâncias, consti- 
tuir-se em transformadas invertíveis que se caracterizam 
justamente em não destruírem informação. 

Este exemplo foi evocado com o fim de mostrar 
como se pode passar de um processo relacional alta- 
mente filtrante e portanto destructivo de informação 
para um processo relacional onde a informação é con- 
servada. 


c) As distribuições foram mencionadas com dois 
fins: 


— Para mostrar uma operação reflexa da da fun- 
cional, 


10 


Com efeito a funcional tem como 1.º mem- 
bro uma relação (função ou distribuição) e por 
2.º membro um real. 

A distribuição de uma variável aleatória 
(que é uma relação), é usada, em controle, como 
2.º membro de uma relação onde o 1.º membro 
é um real (variável determinando potente ou ex- 
plicativa). 


— Para exemplificar como é possível acrescentar 
mais informação a respeito de uma relação cujo 
grafo ocupa um domínio muito vasto para ter in- 
teresse. Foi mostrado que sendo possível proba- 
bilizar o conjunto contradomínio, informação adi- 
cional era susceptível de ser extraída. 

Daí a importância fundamental da def 1'ção 
de medidas sobre espaços mensuráveis quando 
tal é possível. 


3.2 — RELAÇÕES QUE NÃO DESTRÓIEM 
INFORMAÇÃO 


a) Alguns exemplos introdutórios: 


— À corrente eléctrica do par termo-eléctrico tem 
de ser transformada na posição de um ponteiro 
no galvanómetro. 

É essencial que entre a corrente e a posição 
do ponteiro exista uma relação (1 — 1). 


— Às coordenadas de um ponto referidas a duas 
bases, têm de estar de tal modo relacionadas 
que se possa obter uma das outras e inversa- 
mente. 


— Quando numa transformada de Laplace se passa 
do espaço (x, t) para (X, s) é essencial que 
duma função x = qy(t) se passe para um outro 
a X = y(s) e depois seja possível realizar a ope- 
ração inversa e regressar a x = qy(t) sem perda 
de informação. 

Como se vê, a par das Relações que des- 
tróiem informação, tratadas em (3.1), é essencial 
que haja Relações que conservam essa infor- 
mação. 


b) Tem interesse apresentar uma lista de cuidados 
a ter na manipulação da informação de modo a evitar 
que esta se perca, no todo ou em parte. 


1) Podemos concretizar a apresentação do problema 
com o exemplo seguinta: 

O Ponto P num espaço E representa a imagem em 
E de uma certa informação. 

No entanto, por um motivo qualquer é necessario 
ou mais conveniente representar essa imagem no Espaço 
e e o Ponto Q está em correspondência com P. 

Para o efeito usou-se uma relação f tal que 


dom. fc E 
cont. fee 
Às perguntas são: 


— que espécie de espaço e pode servir sabendo 
que a informação original está em E? 
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— que espécie de relação f deve ser usada? 

Se mais nada for dito, a única forma de garantir que 
não haverá, em circunstância alguma, perda de informa- 
ção, é ser-se exremamente exigente: 


— Quanto ao espaço & 
— À dimensão de « e de E tem de ser a mesma. 
— So E tem métrica, e também a tem de ter. 


— Se é possível definir distâncias ou medidas em E 
igualmente terá de ser viável em e. 
isto é, as estruturas algébricas e topológicas de E 
e e têm de ser isomórficas. 


— Quanto a f 


—. f tem de ser uma relação do tipo 1 — 1. 


” “ s pe es 

— Na maioria dos casos, terá fe a sua inversa f - 

de serem contínuas, ou derivadas contínuas até 
qualquer ordem. 


— Em muitos casos será necessária impor que f e 
f * sejam limitados. 


— 6tC., etc. 


Como se vê as condições são tão restritivas cua 
provavelmente, na resolução de problemas reais, podsrão 
concerteza levantar-se algumas destas restrições. 


c) O levantamento de restrições é um problema 
prático e por isso aqui só se poderão apresentar-se al- 
guns casos típicos e ilustrativos. 


— Uma carta de mercator é plana E* é no entanto 
o espiço que procura representar é um x'. À in- 
formação que se porde, os movimentos do navio 
na vertical, não afectando a resolução da classe 
de problemas de navegação marítima para a 
qual as cartas de mercado: foram concebidas, 
torna aceitável este modo de proceder. 


— Um grafo plano das ruas duma cidade onde não 
há passagens sobrepostas é espaço aceitável para 
resolver os problemas de tráfico. 


— Na verdade nestas operações está-se perdendo 
alguma informação, porém essa informação não 
é necessária ou relevante para resolver o Prob- 
lema proposto. 
Assim, a informação inútil para a resolução do 
problema é desprezada. 


4. PROBLEMÁTICA 
4.1 —- ALGUNS PRESSUPOSTOS INTRODUTÓRIOS 

a) Imagem Atributo 

Admite-se um universo de Seres e outro de atri- 
butos. 

Supõe-se que todo o ser tem uma imagem no uni- 


verso dos atributos. 
A imagem «completa» é sempre excessiva em con- 
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teúdo de informação para resolver os problemas do ho- 
mem. 

Recorre-se por isso à imagem «Parcial» adequada 
ao Problema. 

A ciência tem por missão tipificar problemas a es- 
tudar para cada tipo o espaço Imagem adequado. 


b) Restrições 


Os seres interaccionam com outros seres o que per- 
mito ao homem «comandar» ou pelo menos «prevor» 3 
evolução e se o espaço imagem estiver bem construído, 
a imagem evolui à medida que o processo de interacção 
se verifica. 

Uma vez que é possível intervir, podem escolher-se 
os estados permissíveis ao ser (sistema) e assim admi- 
timos a existência de restrições como processo normal 
de intervenção, 

O próprio ser pode ter restrições internas que jus- 
tamente o caracterizam. 


c) Escolha 


Se as restrições internas dos seres e as externas im- 
postas determinarem univocamente um só resultado, o 
problema pode dar-se por terminado. Porém as restrições 
não determinam, em geral, univocamente a solução. 

Neste caso existe um problema de escolha e é es- 
sencial dispor de um critério de escolha. 

O que é escolher? Escolher é muito próximo do 
conceito de ordenar. 

Os métodos mais usados consistem em empregar 
uma relação que tenha por contradomínio o «arquétipo» 
do conjunto ordenado o Corpo dos Reais. 

As vezes uma relação somente, não chega para 9r- 
denar perfeitamente e haverá ainda que escolher ente 
vários elementos ou relações que são equivalentes se- 
gundo o critério de escolha usado. 

Pode repetir-se a operação usando outro critério de 
escolha e assim sucessivamente. 


d) Quem escolhe 


Se existir uma única vontade na definição das res- 
trições e dos critérios de escolha temos o problema mo- 
novolutivo; se intervirem mais do que uma vontade temos 
os problemas plurivolutivos. 

As teorias dos jogos, procuram resolver problemas 
plurivolutivos. 

Podemos dizer que a natureza também exerce a sua 
vontada e os acontecimentos e evoluções que se verifi- 
cam sem intervenção do homem cabem nesta categoria. 

O critério de estremar a entropia em termodinã- 
mica, minimisar a funcional Acção em mecânica clássica, 
são exemplos conhecidos. 

Os critérios não humanos podem igualmente forma- 
lizar-se como acabamos dever no parágrafo anterior e 
então podemos considerar que todo o estado de um ser 
ou a sua evolução resulta do exercício de estremar uma 
funcional da natureza. 

Depois podem ser acrescentadas, nomeadamente, 
as vontades humanas — as quais podem igualmente des- 
crever-se por funcionais extremandas e fronteiras. 


e) Estas funcionais e restrições sobretudo humanas, 
são em geral variáveis e mutantes. 
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4.2 — TIPOS DE PROBLEMAS feridas exercitam-se escolhendo e restringindo. 
As hipóteses de intervenção consideradas no quadro 


Já vimos atrás, que a natureza e os homens podem foram: 
intervir impondo restrições e regras de escolha. — da natúreza 
Normalmente a intervenção exprime-se simultanea- 


mento escolhendo e restringindo, há porém numerosos 
exemplos onde se exerce ou só escolhendo ou só res- 


— de 1 humano 


-— da mais de 1 humano. 


tringindo. O Quadro é suficientemente claro para dispensar 
No quadro junto, admite-se que as intervenções re- comentários. 
Restrições Restrições Restrições 
E) “ e 
Funcionais Funcionais Funcionais 


Naturais 1 Humano > 1 Humano 


Previsões do tempo, do movimento dos astros 
1 0 (0) reacções químicas,etc. 
Funcionais típicas: entropia, Acção, etc. 


Problemas de Comandos, Programa Linear e 
1 1 0) paramétrico, Pert, etc. 


-e Nr ——"|—r -— «a 4 ——— ————— O e 


Jcgos na natureza ou com a natureza. 


1 0 1 

Jogos na natureza na posição de 1 versus 
1 1 1 vários (os mais complexos). 

Caos e experimentação de estruturas sobre 
0 0 (6) o Caos. 

Paciência — jogo para um só jogador que 
(6) 1 0 cria uma regra e a cumpre. 

Jogos entre vários jogadores: 
0 0 1 jogo de cartas, xadres, etc. 

Jogos na posição de 1 contra vários. 
0 1 1 
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DEIXE-NOS RESOLVER O SEU 
PROBLEMA DE FUNDAÇÕES 


Em muitos locais onde vão ser edificadas constru- 
ções sobre terrenos soltos e de fraca consistência, O 
recusrso às técnicas de VIBROFLUTUAÇÃO E 
VIBROCOMPACTAÇÃO é a SOLUÇÃO. 

São técnicas rápidas, que permitem, pela adição e 
vibração em profundidade de areia e/ou brita no 
terreno, a obtenção de uma sólida infraestrutura 
para a construção. 

Os métodos Keller de Vibroflutuação e Vibro- 
compactação conferem ao terreno elevada capaci- 
dade de carga em muitos locais considerados 
antieconómicos para a construção. 


COMPACTE E CONSTRUA 


) 


—- 
— = 


—. 


Assim, antes de tomar uma decisão, escreva ou 
telefone para JOHANN KELLER. O nosso parecer 
pode fazer-vos poupar dinheiro e possibilitar a reali- 
zação dos vossos projectos. 


S. R.L. 


CAPTAÇÕES DE ÁGUA - REBAIXAMENTOS 
AQUÍFEROS - SONDAGENS - CONSOLIDAÇÃO 
DE SOLOS - FUNDAÇÕES ESPECIAIS 


Av. Costa Pinto, 8, 1.º — Cascais 
Tolets. 28 01 46 / 286 44 21 
Telex 16122 
Telegramas WasserK eller 
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PORTO. LERIES 


construção da 


Este empreendimento, adjudicado já depois do 25 de Abril, foi considerado de grande utilidade para a economia do país. 
A natureza dos terrenos obrigou a utilizar duas técnicas distintas para a mesma obra. Assim, num total de 1 200 metros 
de cais acostável 650 serão fundados sobre um tabuleiro apoiado em vigas e estacas de betão armado de grande diâmetro 
(D 1,10 m) e os restantes 580 metros serão constituídos por blocos de betão. 


uma empresa a responder E CONSTRUC OES 


TÉCNICAS, S.A.R.L. 


AV. 24 JULHO, 24 . TEL. 690172 - 691302 - 690163 - LISBOA - 2 


Os necessidades do Pais 
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C. D. U. 532.543.011.1 


Utilização de modelos aerodinâmicos no 


estudo de problemas fluviais 


RESUMO 


A utilização de modelos aerodinâmicos para o es- 
tudo de escoamentos com superfície livre é uma técnica 
que parece ser digna de interesse, por permitir estudar 
certos problemas fluviais em modelos de muito menores 
dimensões do que as dos modelos hidráulicos tradicio- 
nais proporcionando, assim, grandes economias. Por isso, 
foi decidido empreender no LNEC, o estudo dos aspectos 
teóricos e práticos do estabelecimento e exploração dos 
modelos aerodinâmicos, com vista à avaliação das suas 
potencialidades e, eventualmente, à preparação da sua fu- 
tura utilização de rotina. 

Neste artigo faz-se um breve resumo histórico dos 
modelos aerodinâmicos, refere-se o equipamento, mate- 
riais e técnicas de construção a que se recorre, apresen- 
tam-se os critérios de semelhança mais significativos, 
dá-se ideia dos condicionamentos e aplicabilidade desses 
modelos, descrevem-se as operações de taragem e explo- 
ração e comparam-se as vantagens e inconvenientes 
relativos dos modelos aerodinâmicos e dos modelos hi- 
dráulicos tradicionais. Por último relata-se a actividade 
desenvolvida e em curso no LNEC no âmbito dos mo- 
delos aerodinâmicos. 


1 — INTRODUÇÃO 


Como se sabe a utilização de modelos hidráulicos 
para o estudo de problemas fluviais constitui muitas ve- 
zes uma ferramenta da maior utilidade. As dimensões 
destes modelos são condicionadas por um lado pela ex- 
tensão do protótipo que se pretende reproduzir, a qual 
está condicionada, por sua vez, pela natureza do prob- 
lgma a estudar, e por outro lado pelas escalas geométri- 
cas do modelo que são fixadas por forma a satisfazer 
a determinados critérios de semelhança. Na prática os 
modelos hidráulicos fluviais têm muitas vezes dimensões 
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SYNOPSIS 


The use of air-flow models for studying open-chan- 
nel flows seems an interesting technique, because it per- 
mits the study of some river problems with smaller mo- 
dels than the traditional hydraulic ones, thus affording 
great economy. So, it was decided to start at the «Labora- 
tório Nacional de Engenharia Civily the study of theore- 
tical and practical aspects of establishing and using air- 
flow models, to evaluate their potentialities and even- 
tually, to prepare their future use. 

The paper makes a brief historical summary of air- 
flow models, refers to equipment, material and construc- 
tion techniques, shows the most important similar cri- 
teria, gives an idea of the restrictions and applicability 
oi air-flow models, describes the problems of its calibra- 
tion an use and compares the advantages and disadvanta- 
ges of air-flow models with those of traditional hydrau- 
lic ones. Finally it reports on the activity which has been 
carried cut end that is in progress at the LNEC in that 
field. 


apreciáveis que implicam encargos consideráveis e criam 
certas dificuldades na exploração dos modelos. 

Com o objectivo de conseguir modelos com menores 
dimensões foi estudada uma técnica de simulação de 
escoamentos fluviais com superfície livre recorrendo a 
modelos com escoamento de ar em pressão, nos quais 
a superfície livre é simulada por uma placa rígida plana. 
Estes modelos designam-se por modelos aerodinâmicos. 

Na realidade, por razões que adiante se discutirão 
em pormenor, a utilização de ar como fluido no modelo 
o o facto de o escoamento ser em pressão permite que 
os modelos aerodinâmicos tenham dimensões muito me- 
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nores do que os modelos hidráulicos e, consequente- 
mente, sejam mais económicos e em certos aspectos 
mais fáceis de explorar. 

Tanto quanto se sabe, os modelos aerodinâmicos 
foram utilizados pela primeira vez em 1950 no Instituto 
de Investigação de Engenharia Hidráulica de Leninegrado 
na URSS (AVERKIEV 1952) ('). A utilização destes mo- 
delos foi sendo progressivamente aperfeiçoada por di- 
versos investigadores dos países do leste da Europa, no- 
meadamente na própria URSS (AVERKIEV 1954, GILYA- 
ROV 1959, LYATKHER e PRUDOVSKII 1959), na Ro- 
ménia (KAHAN e HÂNCU 1957 e 1958, HÂNCU 1959 
e 1967, HÂNCU e FINKELSTEIN 1963), na Hungria 
(IVICSICS 1966 e 1967), na Checoslováquia (SUMBAL 
1961, 1963, 1964-a e 1964-b, KOMORA e SUMBAL 
1965, KOMORA 1967) e na Polónia (JUNIEWICZ e AU- 
LICH 1958, WALCZAK 1959, JACENKOW e DEMBSKI 
1962). 

Fora destes países apenas se tem conhecimento de 
que tenha utilizado estas técnicas para o estudo de prob- 
lemas fluviais o Laboratório de Hidráulica de Karlsruhe 
na República Federal da Alemanha (BAW 1970) (7). 

Recentemente o Laboratório Nacional de Engenharia 
Civil e o «Laboratoire National d'Hydraulique» de Chatou, 
convictos das possibilidades e do interesse prático dos 
modelos aerodinâmicos, resolveram estabelecer um pro- 
grema de colaboração e empreender em conjunto um es- 
tudo de base com vista a uma futura utilização de rotina 
daqueles modelos. 


2 — CARACTERISTICAS GERAIS DOS MODELOS AERO- 
DINÂMICOS 


Já se referiu que, nos modelos aerodinâmicos, a su- 
perfície livre do escoamento do protótipo é materializada 
no modelo por uma placa rígida plana que deve ser, 
transparente e lisa. Transparente para permitir a visibi- 
lidade e a visualização do escoamento recorrendo a téc- 
nicas fotográficas adequadas. Lisa para assegurar que as 
perdas de carga sejam tão pequenas quanto possível. 

Os modelos aerodinâmicos funcionam normalmente 
em sucção, assegurada por um ventilador colocado na 
extremidade de jusante e têm a extremidade de mon- 
tanto aberta para a atmosfera (Fig. 1). O ar entra por 


“ 
1 I —-—— e 
“POSITIVO DE ; 

| 


) 


esta extremidade e percorre o modelo devido à depressão 
nele criada. Por vezes, colocam-se redes nas extremida- 
des do modelo, não só com a função de evitar a entrada 
de poeiras, mas principalmente para assegurar o amorte- 
cimento das flutuações do escoamento, Antes de atingir 
o ventilador, o ar passa por um dispositivo de medição 
de caudal. 

No caso de modelos com fundo móvel haverá neces- 
sidade de, a montante do modelo, instalar um dispositivo 
de alimentação de caudal sólido e, imediatemente a mon- 
tanto do ventilador, colocar um dispositivo de recolha 
das particulas transportadas através do modelo. 

As condutas e transições podem ser de qualque: 
material, embora a utilização de materiais plásticos apre- 
sente a vantagem de provocar perdas de carga reduzidas. 

A regulação do caudal de ar no modelo pode ser 
feita por meio de um regulador da velocidade do motor 
de um variador de velocidade na transição entre o motor 
e o ventilador ou à custa de uma comporta colocada 
junto ao ventilador o que constitui a solução mais econó- 
mica. 

Na medição do caudal de ar escoado nos modelos 
podem-se utilizar um aparelho deprimogéneo (diafragma 
ou venturi), um tubo de Pitot fixo ou um simples conta- 
dor de gás. 

O ventilador adequado para estas instalações é do 
tipo radial corrente no mercado. É aceitável a utilização 
de um motor eléctrico com velocidade constante, po- 
dendo fazer-se a sua montagem acopulada ao eixo do 
ventilador ou recorrer à utilização de transmissões por 
correias. Para escolher o ventilador é necessário o conhe- 
cimento aproximado das perdas de carga do escoamento 
ao longo de todo o sistema para a gama de caudais de 
funcionamento do modelo. A decisão a tomar deverá ter 
em consideração estes elementos e as curvas caracte- 
rísticas de cada ventilador. Por vezes, quando as perdas 
du carga são muito elevadas, torna-se necessário utilizar 
um compressor em substituição do ventilador. 

Para se conseguir uma precisão suficiente na repro- 
dução da geometria do protótipo é necessário ter muito 
cuidado na construção dos modelos aerodinâmicos em 
consequência das reduzidas escalas utilizadas. Assim 
os materiais, os instrumentos de moldagem e as técnicas 
de execução devem ser escolhidas por forma a permitir 
a precisão desejada. 


] |) 
1 DISPOSITIVO DE t 
+ ALIMENTAÇÃO DE MODELO ANTA 
: Í MEDIÇÃO DE CAUDA.S ] à 
| CAUDAL SÓLIDO | DE PARTICULAS | 
] ] 


equipamento para qualquer tipo de modelo 


.——— - == equipamento específico dos modelos com fundo movel 


Fig. 1 ——- Esquema da instalação de um modelo aerodinâmico 


(1) — As referências constam da bibliografia que se apresenta no final. 


(2) — Também se tem notícia de um caso de utilização de modelos aerodinâmicos para o estudo de problemas hidráulicos no Insti- 
tuto de Investigação Hidráulica da Universidade de lowa nos Estados Unidos da América o qual, contudo, não se refere a problemas 


fluviais (ROUSE 1947). 
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Os modelos com fundo fixo devem ser construídos 
com um material facilmente moldável e que assegure 
a estanquidade. Usa-se com frequência a plasticina, o 
gesso, a parafina e a estearina. Nos modelos com fundo 
móvel têm sido utilizados para reproduzir o fundo mate- 
riais muito leves que possam ser transportados pelo ar, 
sendo corrente usar areia fina, carvão ou farinha. Parece 
provável, contudo, que, entre a vasta gama de materiais 
sintéticos de que se dispõe hoje em dia, se possam en- 
contrar alguns particularmente recomendáveis para este 
fim. 

A placa que cobre superiormente os modelos deve 
ser de vidro ou de plástico transparente, podendo ter ori- 
fícios em diversas secções transversais para a introdu- 
ção dos aparelhos de medida. O vidro origina menores 
flechas que o plástico mas, por outro lado, é mais que- 
bradiço, sendo difícil a abertura dos orifícios. Para evitar 
flechas inconvenientes da placa, pode-se apoiar esta em 
pequenos apoios situados no interior do modelo ou ligar 
a placa a uma estrutura exterior. 

É essencial garantir a completa estanquidade de 
toda a instalação, Assim, todas as juntas e todos os furos 
existentes na placa quando não estejam a ser utilizados 
devem ser perfeitamente obturados. Este problema é de 
certo modo facilitado por o modelo funcionar em sub- 
pressão. Normalmente para obturar os furos utilizam-se 
pequenas tampas de borracha, plástico, cortiça ou outro 
material e, nas juntas, usam-se tiras de borracha com pe- 
quena espessura fortemente comprimidas entre as peças 
a unir que são, em geral, aparafusadas entre si. Para 
garantir maior eficiência à vedação pode recorrer-se a 
fita adesiva. 

Os modelos aerodinâmicos são tão pequenos que 
é possível, em geral, utilizá-los sobre mesas. Os modelos 
de fundo fixo podem, inclusivamente, funcionar pendu- 
rados numa parede, pois nos escoamentos em pressão é 
desprezável a acção da gravidade. Deste modo o espaço 
ocupado pelos modelos aerodinâmicos é ainda mais re- 
duzido. 


3 — CONDICIONAMENTOS DOS MODELOS 
AERODINÂMICOS 


As escalas lineares dos modelos hidráulicos para 
o estudo de problemas de escoamentos com superfície 
livre são normalmente fixadas com base nos dois condi- 
cionamentos seguintes: 


— necessidade de o escoamento no modelo se pro- 
cessar em regime turbulento rugoso. 


— necessidade de respeitar o critério de seme- 
lhança de Froude. 

O primeiro condicionamento provém da necessidade 
de assegurar que no modelo, tal como sucede no pro- 
tótipo, as perdas de carga sejam proporcionais ao qua- 
drado da velocidade. A garantia de o escoamento ser 


turbulento rugoso é assegurada pela condição (SCHLI- 
CHTING 1960): 


V,H 
Re =——— 70 
às o? (3.1) 
onde: 
= V —t  — velocidade de atrito no fundo 
Ê 


H — profundidade média do escoamento 


r — viscosidade cinemática do fluido utilizado 
no modelo 


ct — tensão tangencial média no fundo 


p — massa específica do fluido 


O segundo condicionamento resulta de, na maioria 
dos problemas fluviais haver preponderância das forças 
da gravidade sobre as forças de viscosidade e tem como 
consequência a igualdade entre a escala geométrica à y 
e o quadrado da escala das velocidades Ay (3) 


Ap=A “v (3.2) 

As dimensões mínimas dos modelos hidráulicos são 
deste modo condicionadas pela necessidade de respeitar 
os dois condicionamentos referidos. 

Uma forma de reduzir as dimensões de um modelo 
hidráulico seria utilizar no modelo um fluído com menor 
viscosidade cinemática do que a da água. No entanto 
esto procedimento, em face da gama relativamente redu- 
zida das viscosidades dos fluídos disponíveis nunca per- 
mitíria reduzir muito as dimensões do modelo e, além 
disso, não é económico. 

A utilização de modelos aerodinâmicos permite con- 
seguir o regime turbulento rugoso com escalas geomé- 
tricas muito menores do que as necessárias para se con- 
seguir o mesmo regime dos modelos hidráulicos, porque 
nos modelos em pressão é possível obter velocidades 
muito mais elevadas do que nos modelos com superfi- 
cio livre e, além disso, não se está condicionado pela se- 
melhança de Froude. Deste modo, a utilização de veloci- 
dades elevadas não obriga a construir os modelos com 
grandes dimensões para respeitar a condição imposta 
pela equação (3.2). Note-se que as velocidades que se 
utilizam nos modelos aerodinâmicos são suficientemente 
elevadas para que a condição (3.1) possa ser verificada, 
apesar de a viscosidade cinemática do ar ser, para as 
temperaturas correntes de ensaio, maior do que a da 
água, tal como se pode ver na Fig. 2 onde se mostra a 
variação da viscosidade cinemática do ar e da água com 
a temperatura. 

A velocidade no modelo deve ter um limite supe- 
rior, determinado pela compressibilidade do ar e por 
condicionamentos inerentes à exploração do modelo. 


(1) Representa-se a escala de uma grandeza X por Ax como sendo a razão entre os valores que a grandeza X toma no modelo 


Xm 
Xp 


(Xm) e no protótipo (X,): Ay = 
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Osho 


VISCOSIDADE CINEMATICA v (m2/5) 


AJA tetead= ES 
-40 - 0 o 4 re 190 
TEMPERATURA & [ºC] 


Fig. 2 — Influência da temperatura na viscosidade cinemática do 
ar e da água 

Para se determinar o valor máximo da velocidade 
pode recorrer-so à análise da influência da velocidade 
nas diferentes características do ar. Na Fig. 3 está re- 
presentada a variação da carga dinâmica com a veloci- 
dade do ar considerado compressível e incompressível, 
o que permite concluir que a partir das velocidades da 
ordem de 100 m/s a compressibilidade do ar começa 
a introduzir erros significativos. Alguns autores apontam 
para valor limite superior da velocidade do ar valores 
ente 50 e 80 m/s. Por razões ligadas com dificuldades 
da exploração dos modelos e em particular com a garan- 
tia de condições de estanquidade, os valores máximos 
da velocidade utilizados na prática são geralmente bas- 
tante inferiores aos condicionados pela compressibilidade 
do ar, não ultrapassando em regra valores da ordem 
de 20 m/s. É conveniente notar que o valor limite supe- 
rior não deve ser atingido em nenhum ponto do modelo, 


CARGA DINÂMICA 


350 
YELOCIDADE Vims) 


Fig. 3 — Influência da compressibilidade do ar na relação entre 
a carga dinâmica e a velocidade média 
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o que obriga a que ele seja respeitado na secção de me- 
nor área. 


4 — CRITÉRIOS DE SEMELHANÇA UTILIZADOS NOS 
MODELOS AERODINÂMICOS 


Apesar de haver a preocupação de que a placa que 
simula a superfície livre nos modelos aerodinâmicos seja 
tão lisa quanto possível, geram-se junto à placa esforços 
tangenciais maiores do que os correspondentes no protó- 
tipo à superficie de contacto ar-água, os quais são, em 
geral, desprezáveis (Fig. 4a). 


EPP SPP 


Fig. 4 — Distribuições verticais de velocidade: a) protótpo; b) 
modelo sem sobreelevação da placa (semelhança global); c) mo- 
delo com sobreelevação da placa (semelhança local) 


Por consequência, a estrutura cinemática do escoa- 
mento é alterada e os diagramas de distribuição vertical 
de velocidades não se correspondem (Fig. 4 b). 

Nos modelos aerodinâmicos adoptam-se critérios 
de semelhança definidos de acordo com o tipo de prob- 
lemas a estudar. Quando se pretende apenas conhecer 
o andamento da superfície livre ou o valor médio de gran- 
dezas características do escoamento aplicam-se critérios 
do semelhança global. Quando além da superfície livre 
e dos valores médios das grandezas características do 
escoamento se pretende conhecer os seus valores locais 
Be, em particular, caracterizar a sua estrutura cinemática 
(distribuição de velocidades, campo de linhas corrente, 
etc.) aplicam-se critérios de semelhança local. 

No primeiro tipo de problemas respeita-se a seme- 
lhança geométrica, sendo a placa superior colocada no 
modelo na posição corespondente à superfície livre do 
escoamento do protótipo, e despreza-se a semelhança 
entre os diagramas das velocidades (Fig. 4 b). 

Nos problemas do segundo tipo sacrifica-se a seme- 
lhança geométrica, sendo a placa elevada para uma po- 
sição tal que permita que a superfície das velocidades 
máximas do modelo fique a corresponder, por seme- 
lhança geométrica, à superfício livre do protótipo, Ga- 
rante-se assim a analogia entre a distribuição de veloci- 
dades do protótipo e a distribuição de velocidades na 
parte inferior do modelo (Fig. 4 c). 

Analisa-so em seguida como se costumam estabe- 
lecaer os critérios de semelhança global utilizados nos 
modelos aerodinâmicos. A partir da equação de energia 
aplicada ao escoamento no protótipo e ao escoamento 
no modelo é possível chegar às seguintes igualdades: 
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SP m 
q Va SE, 
Az, Vo AE, (4.1) 
onde: 
á ty diferença da cota da superfície livre entre duas 
secções do protótipo 
à Pp, — diferença de pressão estática entre as duas sec- 


ções correspondentes do modelo 
V — velocidade média do escoamento 


à E — perda de energia do escoamento por unidade de 
peso entre duas secções de referência 
yY — peso específico do ar 
representando os índices m e p as grandezas que se re- 
ferem ao modelo e ao protótipo respectivamente. 
A partir da primeira igualdade de (4.1) e fazendo 
y = Bg chega-se ao seguinte parâmetro adimensional: 


95Zz, SP, 
a “Wu 42) 
V p p m (4. 
onue: 
Eu — número de Euler 
: — massa específica do ar 
g — aceleração da gravidade 
Uma vez que 
IV: 
dA E ——— 
8 g R (4.3) 


onde: 
f —- coeficiente de resistência ao escoamento de 
Darcy-Weisbach 


| — distância entre as duas secções consideradas 
V — velocidade média do escoamento 
R — raio hidráulico 


da ser «ua igualdade de (4.1) resulta: 


2 m Dm R,p = 
fp Ip Rm (4.4) 
Lu 
ke À 
É o 
ÀpR (4.5) 


Ora, só um dos critérios (4.2) e (4.5) é determi- 
nante. HÂNCU 1967 considerou como determinante o 
segundo o que corresponde a aceitar que a semelhança 
entre os dois escoamentos fica assegurada se se verifi- 
car (4.5) sendo o parâmetro adimensional (4.2) função 


51 
do parametro: , OU Seja: 


fl 
Eu=q (—) 
R 
No caso particular em que se considerem canais su- 
ficientemente largos a relação de escalas 14,5) pode es- 
crever-se: o 


(4.6) 


= 1 (4.7) 
An 


onde h é a profundidade do escoamento. 
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So além de se considerarem canais suficientemente 
largos o modelo não for distorcido ter-se-á ainda 

Ap =1 (4.8) 
o que implica que neste caso o factor de resistência no 
protótipo seja igual ao factor de resistência no modelo. 

Atendendo a que o escoamento é uniforme e con- 
siderando a fórmula de Strickler 


A E 

V=kK R35il2=KçR 25 ( ) 1/2 ai 
onde: 

K. — coeficiente de rugosidade global de Strickler 

Do. perda de carga por unidade de peso e de per- 
curso 

vem 
) E o 
= —— 
a (4.10) 


Tendo em conta que (4.5) pode ser expresso por: 


VE Ap (4.11) 
que se obtém directamente de (4.1), e combinando 
(4.10) e (4.11) obtém-se: 


(4.12) 


Os critérios de semelhança local podem ser também 
estabelecidos de forma análoga à anterior, bastando ter 
presente que se aplicam apenas à parte inferior do mo- 
delo. 

No entanto GILYAROV 1959 e MAKKAVEEV 1959 
deduziram critérios de semelhança local a partir das 
seguintes igualdades: 


p lim (4.13) 
válidos para canais suficientemente largos em que se 
possa considerar o raio hidráulico igual à profundidade 
é sendo: 


Hp o Ma inclinação longitudinal respectivamer:oe 


da superfície livre do protótipo e da 
placa superior do modelo 
Hp é "em — inclinação transversal respectivamente 


da superfície livre do protótivo e da 
placa superior do modelo 


GILYAROV 1959 conclui que, no caso da elevação 
da placa obedecer à relação 


, Ei ha (4.14) 
onde: 
hm — profundidade total do modelo 
h;m — profundidade da parte inferior do mo leto 
a escala de velocidades é dada por: 
o dim o (415) 
onde: 


Yi 


Ah; — escala das profundidades considerando no 
modelo apenas a profundidade da sua parte 
inferior. 


Num caso mais geral em que se eleve a placa sem 
obedecer a (4.14) e também segundo GILYAROV 1959 
a escala de velocidades poderá ser calculada pela re- 
Izção: 


him o 
NG A E ) (4.16) 


onde: 


C —- coeficiente de Chézy. 


5 —. TARAGEM E EXPLORAÇÃO DOS MODELOS 
AERODINÂMICOS 


5.1 —- GENERALIDADES 


À semelhança do que acontece quando se utilizam 
modelos hidráulicos, também a reprodução correcta em 
modelos aerodinâmicos dos fenómenos que ocorrem no 
protótipo não é automaticamente garantida pela constru- 
ção e exploração dos modelos resultantes da aplicação 
dos critários de semelhança a que na prática se recorre. 
É necessário fazer a taragem dos modelos, isto é, simu- 
lar nos modelos, construídos de acordo com os critérios 
de semelhança adoptados, fenómenos que já tenham 
ocorrido no protótipo. 

Pode verificar-se assim, se o modelo reproduz ou 
não correctamente o protótipo no que se refere aos fenó- 
menos a estudar e, caso não reproduza, torna-se neces- 
sário modificar, por tentativas, a rugosidade do modelo 
atá se obter a rugosidade correcta. Esta taragem, efec- 
tuada por aproximações sucessivas, é do mesmo tipo da 
que se realiza nos modelos hidráulicos. Depois de a tara- 
gem estar concluída pode passar-se à exploração, en- 
saiando as obras que se pretendem estudar. 

Faz-se notar que tanto na taragem como na explora- 
ção dos modelos se deve medir a temperatura do ar, 
a pressão e a humidade atmosféricas para fazer as neces- 
sárias correcções das grandezas afectadas, nomeada- 
mente do peso específico do ar. 

Nas duas secções seguintes descrevem-se a taragem 
e a exploração de modelos aerodinâmicos de fundo fixo 
quando interessa apenas o estudo de grandezas médias, 
isto é, seguindo critérios de semelhança global. Não se 
faz referência à taragem e à exploração dos modelos com 
fundo fixo quando são seguidos os critérios de seme- 
lhança local, nem dos modelos com fundo móvel, por 
não se conhecerem ainda critérios suficientemente ex- 
plícitos e inequívocos. 


5.2 — TARAGEM DOS MODELOS AERODINÂMICOS 


meira equação da expessão (4.1) que se pode escrever: 

Analisa-so em seguida o procedimento normal para 
realizar a taragem de um modelo aerodinâmico concebido 
poi forma a satisfazer a critérios de semelhança global. 
Aceita-so o critério de semelhança traduzido pela pri- 
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ap 


m 
; 
AZ, 4 (5.1) 


onde (Figs. 5 e 6): 


| 
ais 


—. diferença das cotas da su- 
perfície livre no protótipo, 
entre as secções 1 6 2. 


az, =2Zp op 


SP =Pia "Pam —. diferença entre valores da 
pressão no modelo entre as 
secções correspondentes, 1 
e 2. 


Y —peso específico do fluído 
no modelo. 


Uma vez que são conhecidas as características 
geométricas do protótipo e as cotas da superfície livre 
nas secções 1 e 2 para cada caudal, é possível construir 
o modelo obedecendo a determinadas escalas geométri- 
cas e colocar a placa superior na posição correspondente 
à superfície livre. 

A taragem faz-se, tal como nos modelos hidráulicos, 
à custa da variação da rugosidade do fundo. A primeira 
rugosidade a ensaiar vode ser arbitrada, mas, no caso 
particular de canais muito largos e desde que se conheça 
a rugosidade do protótipo é possivel estimar um primeiro 
valor da rugosidade do modelo recorrendo à equação 
(4.12). 

Depois de intoduzida no modelo a rugosidade es- 
colhida, faz-se escoar um caudal de ar que apenas tem 
de obedecer aos condicionamentos gerais indicados em 
3. Em geral utilizam-se caudais que correspondam a ve- 
locidades da ordem de 20 m/s pelas razões já referidas 
em 3. 

Nesta altura é possível determinar pela equação 
(5.1) o valor de Ap, . pois y e à Z, são conhecidos e 
ày pode ser calculado a partir dos valores da velocidade 
conhecida no protótipo e medida no modelo. Como, por 
outro lado, se pode medir no modelo esse mesmo valor 


AP, apenas há que comparar o valor de Ap, calculado 


(4 pf ) com o valor de Ar, determinado experimen- 
talmente (4 Pim ). Esta comparação pode conduzir a duas 
situações: 


1) [À Pm — A Pg | > 81 


E e 
2) |4pm —SAP| Sô: 


sendo 5, o erro máximo admissível para a diferença de 
pressão entre as secções 1 e 2 do modelo. 

A situação 1) corresponde a aceitar que a equação 
(5.1) é verificada, considerando-se, portanto, a taragem 
realizada como satisfatória. 

Seo se verificar a situação 2) podem ainda dar-se 
dois casos diferentes: 


2 a) APm dA 
2 b) Apm CAPA: 


No primeiro caso a perda de carga do modelo é de- 
masiado pequena, sendo necessário aumentar a rugosi- 
dade e repetir o ensaio até que se verifique a situação 1). 
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No segundo caso é necessário diminuir a rugosidade até 
que a situação 1) se verifique. 

A taragem é válida para qualquer caudal, havendo 
no entanto vantagem em mantê-lo constante nas diversas 
fases da taragem. Assim evita-se ter de, para cada ten- 


tativa, recalcular o valor à Pm já que deste modo hy se 
mantém constante, 


5.3 — EXPLORAÇÃO DOS MODELOS AERODINÂMICOS 


Para exemplificar como se faz a exploração de um 
modelo aerodinâmico considera-se também o caso de um 
modelo concebido por forma a satisfazer os critérios de 
semelhança global. O modelo representa-se esquematica- 
mente na Fig. 5 b, estando a secção 1 suficientemente 
afastada da zona em que se pretende construir uma obra 
a ensaiar para, por ela, não ser influenciada. Pretende-se 
utilizando o modelo, prever, para um determinado caudal, 
a cota Z'op que terá a superfície livre do protótipo na 
secção 2 depois de construída uma obra, por exemplo 
de correcção fluvial, que está representada esquematica- 
mente na Fig. bc. Para isso recorre-se, tal como na ta- 
ragem, a um método de aproximações sucessivas que 
permite determinar a posição da placa superior que cor- 
responde, de acordo com a escala geométrica adoptada 
para o modelo, à superfície livre do escoamento no pro- 
tótipo após a construção da obra. 

Representa-se então no modelo a obra e coloca-se 
a placa numa posição que a priori se supõe corresponder 
à nova posição da superfície livre no protótipo, tendo no 
entanto o cuidado de manter a cota na secção 1, isto é, 
fazendo Z', | =Zym - À cota da placa na secção 2 desig- 


na-se por Z' m (1) à que corresponde no protótipo Z'o p(1) 

Para um determinado caudal que se faça escoar no 
modelo e que obedeça às condições de funcionamento 
atrás referidas, é possível, atendendo aos valores conhe- 
cidos na secção 1, detarminar Av Esta determinação 
já estará feita anteriormente se se mantiver o caudal uti- 
lizado na taragem. 

De posse destes elementos e medindo a diferença 


AZprzip-r2 


de pressão entre as duas secções de referência, a apli- 
cação da equação (5.1) permite calcular analiticamente 
um valor da cota da superfície livre na secção 2 após a 
obra construída que se designa por Z'op (9) € à que cor- 
responde no modelo uma cota Zám (2) 


Ão comparar z op (1) COM Z'95 (9) dois casos podem 
ocorrer: 


Dz cp(D —2op(2) | 2 8 
2) | z op(1) Top (9) |S ô 


onde 3. é o erro máximo admissível para a cota da su- 
perfície livre na secção 2 do protótipo. 

No caso 1) os valores são suficientemente próximos 
para se concluir que a equação (5.1) é verificada e por- 
tanto a posição em que se encontra a placa representa 
correctamente a superfície livre do escoamento que ocor- 
rerá no protótipo se a obra for realizada. Assim fica 
conhecida a cota Z'op como se pretendia. 

No caso 2) os valores não são suficientemente pró- 
ximos para que a equação (5.1) se considere verificada 
e, portanto, é necessário alterar, por aproximações suces- 
sivas, a posição da placa até se conseguir verificar o 
caso 1) ou, o que é mesmo, satisfazer a equação (5.1). 


6 — COMPARAÇÃO DOS MODELOS AERODINÂMICOS 
COM OS MODELOS HIDRÁULICOS 


As vantagens dos modelos aerodinâmicos relativa- 
mente aos modelos hidráulicos derivam de uma substân- 
cial economia ds material, de tempo e de espaço. Esta 
economia resulta, por um lado, de se utilizarem modelos 
com dimensões consideravelmente inferiores às dos mo- 
delos hidráulicos e, por outro, de se utilizar como fluido 
o ar em vez da água. 

As principais vantagens dos modelos aerodinâmicos 
relativamente aos modelos hidráulicos podem-se enu- 
merar da seguinte forma: 


— menor ocupação de espaço, o que permite que a 
área dos laboratórios destinada à instalação dos modelos 


Pm * Pim” P2m , 


dedo 


Fig. 5 —- Representação esquemática do escoamento em corte longitudinal: a) protótipo (escoamento da água com superfície 
livre); b) modelo (escoamento de ar em pressão); c) idem, com obra implantada no leito 
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seja consideravelmente reduzida e possibilita que se re- 
produzam nos modelos maiores troços dos rios a estudar; 


— menor quantidade de material e mão-de-obra 
utilizados tanto na construção dos modelos como nas mo- 
dificações que se têm de introduzir durante os ensaios; 


— maior possibilidade de recuperar, quando da 
destruição dos modelos, grande parte do material para 
construção de novos modelos; 


— maior economia na instalação, por se dispensa- 
rem circuitos hidráulicos, com reservatórios de montante 
e de jusante, e sistemas de bombagem mais ou menos 
complexos; com efeito nos modelos aerodinâmicos os 
reservatórios são o próprio ar atmosférico e o sistema 
de bombagem é um simples ventilador accionado por 
um motor eléctrico com potência relativamente pequena; 


— maior comodidade e rapidez na exploração dos 
modelos; 

— maior possibilidade de ter uma visão global dos 
fenómenos que ocorrem nas zonas a estudar. 


No caso de modelos com fundo móvel os modelos 
aerodinâmicos permitem ainda, relativamente aos mode- 
los hidráulicos, as seguintes vantagens: 

— menor quantidade de material aluvionar utilizado; 


— maior facilidado na exploração e consequente- 
mente redução da duração dos ensaios; 


— maior simplicidade nas operações de moldagem 
6 levantamento da topografia dos fundos. 


Os inconvenientes dos modelos aerodinâmicos re- 
lativamente aos modelos hidráulicos resultam, na sua 
maior parte, da necessidade de materializar a superfície 
livre do escoamento do protótipo por uma placa rígida 
e da maior delicadeza de que se reveste a construção dos 
modelos e o trabalho experimental. Estes inconvenientes 
dos modelos aerodinâmicos relativamente aos modelos 
hidráulicos podem-se enumerar da forma seguinte: 


— menor correcção na representação do escoa- 
mento do protótipo por existirem, junto à placa que si- 
mula a superfície livre, tensões tangenciais que alteram 
a estrutura cinemática do escoamento e afectam o gra- 
diente de velocidades; 


— impossibilidade de reprodução de escoamentos 
rapidamente variados devido à presença da placa que tem 
uma posição fixa; 

—— necessidade de se utilizarem instrumentos de 
medida mais precisos e sensíveis e exigência de maior 
cuidado e precisão na construção dos modelos. 


Pelo que se disse, vê-se que os inconvenientes mais 
correntemente apontados aos modelos hidráulicos, tais 
como as grandes dimensões, a necessidade de dispor de 
grandes caudais, os elevados custos e os dilatados pra- 
zos de realização de ensaios, são muito satisfatoriamente 
superados nos modelos aerodinâmicos, embora, por ve- 
zes, à custa de uma certa redução de precisão. No en- 
tanto, em grande numero de casos verifica-se que os 
modelos aerodinâmicos dão resultados suficientemente 
precisos e permitem substituir com grande economia de 
custos e de tempo os modelos hidráulicos. Além disso, 
mesmo quando os modelos hidráulicos não possam ser 


completamente dispensados, os modelos aerodinâmicos 
constituem um meio suficientemente seguro para uma 
primeira selecção das hipóteses inicialmente encaradas 
ou permitem a definição das condições de fronteira a uti- 
lizar posteriormente em modelos hidráulicos que repro- 
duzam áreas mais circunscritas representadas a escalas 
maiores, 


7 — APLICABILIDADE DOS MODELOS 
AERODINÂMICOS 


Como se viu em virtude de nos modelos aerodinã- 
micos a superfície livre do escoamento no protótipo ser 
simulada por uma placa rígida, estes modelos são ade- 
quados para o estudo de problemas em que o escoamento 
é permanente, não se aplicando normalmente a casos 
de escoamento variável. Além disso, uma vez que a placa 
além de ser rígida é, em princípio, plana, os modelos 
aerodinâmicos são particularmente úteis para o estudo 
de escoamentos permanentes gradualmente variados, cuja 
superfícies livre se pode considerar praticamente plana 
para as extensões reproduzidas nos modelos. 

Os modelos aerodinâmicos permitem estudar tanto 
escoamentos com fundo fixo como escoamentos com 
fundo móvel. No entanto, no estudo dos escoamentos 
com fundo móvel a natureza dos resultados que é pos- 
sível obter é diferente no caso da fase líquida e no da 
fase sólida. 

Na realidade, a fase líquida é susceptível de melhor 
caracterização, uma vez que se podem obter resultados 
tanto de carácter qualitativo, como, por exemplo, iden- 
tificação de linhas de corrente e turbilhões, como de ca- 
rácter quantitativo, como, por exemplo, medição de velo- 
cidades pontuais, de caudais e de pressões. 

Relativamente à fase sólida os resultados obtidos 
são meramente qualitativos, embora seja suficientemente 
correcta a reprodução de fenómenos morfológicos locais 
ou globais e, portanto, das erosões e assoreamentos, da 
meandrização, etc. 

Embora haja um grande número de problemas que 
podem ser completamente resolvidos pelo recurso aos 
modelos aerodinâmicos, dispensando-se, portanto, a 
construção de modelos hidráulicos, outros existem em 
que os resultados obtidos com base nos modelos aerodi- 
nâmicos não são, só por si, suficientemente correctos, 
devido principalmente à deformação das características 
cinemáticas do escoamento ocasionada pela placa, ha- 
vendo portanto vantagem em recorrer aos modelos hi- 
dráulicos. 

Porém, mesmo neste último tipo de problemas, Jus- 
tifica-se muitas vezes a construção simultânea de um mo- 
delo hidráulico e de um modelo aerodinâmico pela 
grande economia que se consegue. Por um lado, a cons- 
trução de um modelo aerodinâmico que abranja uma re- 
gião mais vasta no protótipo permite obter, em pouco 
tempo, condições de fronteira correctas para um modelo 
hidráulico que reproduza uma zona mais restrita do pro- 
tótipo. Por outro lado, o modelo aerodinâmico pode per- 
mitir ensaiar rapidamente um grande número de alterna- 
tivas para a resolução do problema em estudo, Assim, 
só as alternativas que se revelem mais convenientes 
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serão posteriormente ensaiadas num modelo hidráulico 
em maior escala para se escolher a solução definitiva 
para o problema. 


8 — ACTIVIDADE DO LABORATÓRIO NACIONAL DE 
ENGENHARIA CIVIL RELACIONADA COM O ES- 
TUDO DOS MODELOS AERODINÂMICOS 


Por se considerar que os modelos aerodinâmicos 
podem constituir uma preciosa ferramenta para o estudo 
dos problemas fluviais, muito útil e insuficientemente 
explorada em todas as suas potencialidades, o LNEC 
decidiu promover um estudo sistemático desta nova téc- 
nica com vista a esclarecer certos aspectos do funciona- 
mento dos modelos aerodinâmicos e a possibilitar a sua 
futura utilização de rotina. Este estudo foi lançado em 
colaboração com o «Laboratoire National d'Hydraulique» 
de Chatou que também estava interessado no problema. 

Previu-so que o trabalho a realizar no LNEC fosse 
executado nas quatro fases seguintes: 


1.º fase: — apropriação de conhecimentos 

2.º fase: — aquisição de aparelhagem e construção 
de instalações piloto 

3.* fase: — esclarecimento de diversos problemas 
recorrendo às instalações piloto 


4.º fase: — duplicação de um caso concreto estu- 


dado em modelo hidráulico. J 


o. 


Na 1.º faso procedeu-se a uma exaustiva pesquisa 
bibliográfica e traduziram-se alguns artigos e obras que 
se consideraram mais significativos (visto que grande 
parto dos textos fundamentais está publicada em idiomas 
não acessíveis) tendo-se elaborado um relatório de sín- 
teso que constitui o documento que serve de base para 
o prosseguimento do trabalho. 


Ainda nesta fase foi realizada uma visita de estudo 
extremamente proveitosa ao «Instituto de Investigação 
para a Recuperação de Terras e a Gestão das Águas» 
de Bucareste na Roménia, ao «Instituto de Investigação 
para o Desenvolvimento dos Recursos Hídricos» (VITUKI) 
de Budapeste na Hungria, ao «Instituto de Investigação 
de Hidrodinâmica» (VUV) de Bratislava na Checoslová- 
quia e ao «Instituto Federal para as Obras Hidráulicas» 
(BAW) de Karlsruhe na Alemanha Federal, laboratórios 
e equipamento necessários para a exploração dos mo- 
delos aerodinâmicos e construiram-se instalações piloto. 


Na 2.º fase procedeu-se à aquisição de aparelhagem 
e equipamento necessários para a exploração dos modelos 
aerodinâmicos e construiram-se instalações piloto. 

Entro a aparelhagem de medição adquirida ou cons- 
truída nas oficinas do LNEC inclui-se um conjunto com 
termómetro, higrómetro e barómetro, dois micromanó- 
metros simples de álcool inclináveis, um multimanómetro 
de álcool inclinável, duas sondas para medição de velo- 
cidades (um tubo de Pitot e uma sonda esférica tipo 
esfera de Pitot) e um suporte para sondas permitindo 
medir as rotações e deslocamentos verticais das sondas 
(Fig. 6). 

As instalações piloto compreendem um conjunto 
fixo a jusante, constituído por dispositivo de medição 


Fig. 7 —- Aspecto do troço ds transição de jusante, sistema de 
medição de caudais e ventilador 
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de caudais, ventilador e transições de secção (Fig. 7), 
a montanto do qual é possível adoptar diferentes mo- 
delos. Os dois modelos até agora construídos (Fig. 8) 
têm secção rectangular com altura regulável em virtude 
dao a placa superior ser móvel. Esta nlaca tem uma série 
do orifícios para introdução das sondas de medição dis- 
postos em diversas secções. Os dois modelos diferem 
apenas em planta pois um deles é rectilíneo e o outro 
tem forma de U. 


a 


A , a 
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Fig. 8 —- Aspecto dos mode'os rectilineo e em U 


A 3.º fase, que está em curso, consta de um traba- 
lho experimental com o objectivo de verificar a validade 
dos conceitos teóricos em que se fundamenta a concep- 
ção, exploração e taragem dos modelos aerodinâmicos, 
de ensaiar a aparelhegem de medida e de tomar contacto 
com as dificuldades que possam surgir na exploração dos 
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modelos. Um dos pontos que requer especial atenção 
é o problema da reprodução de correntes secundárias 
nos modelos aerodinâmicos, em particular as resultantes 
dos escoamentos em curvas. Para a realização deste es- 
tudo tem particular interesse a instalação piloto curvilí- 
nea a que se fez referência. 

Paralelamente decorre o estudo teórico, incidindo 
especialmente no estudo dos critérios de semelhança 
a adoptar para simular escoamentos de água com super- 
fícia livro por meio de escoamentos de ar em pressão, 
pois os autores que estudaram o problema não estão 
completamente de acordo quanto aos critérios de seme- 
lhança a utilizar nem os fundamentam devidamente numa 
perspectiva racional. 

Por fim, na 4.º fase, prevê-se a construção de um 
modelo aerodinâmico duplicando um caso concreto que 
esteja a ser estudado em modelo hidráulico, para assim 
se poderem comparar as medições de exploração e os re- 
sultados dos modelos aerodinâmicos com os dos modelos 
hidráulicos. 

Após o termo desta última fase pensa-se poder 
começar a utilizar os modelos aerodinâmicos com carác- 
ter de rotina nos numerosos problemas em que a sua uti- 
lização parece ser vantajosa. 
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Sobre a aplicabilidade dum modelo à previsão 


da poluição atmosférica 


RESUMO 


Este artigo destina-se a pessoas que, não estando 
muito familiarizadas com modelos de previsão de polui- 
ção atmosférica, estejam mais interessadas em conhecer 
as suas potencialidades do que as bases físicas do pro- 
cesso. 

Descreve em termos gerais as utilizações mais cor- 
rentes destes modelos nomeadamente planificação da 
utilização do solo, avaliação de estratégias de controle 
e funcionamento de sistemas de alerta. 

Os modelos Gaussianos, apesar das suas limitações, 
têm permitido obter resultados satisfatórios, são de uti- 
lização fácil e necessitam pouco tempo de computador, 
pelo que lhes é dado especial relevo. Referem-se os 
dados necessários ao funcionamento destes modelos 
assim como as informações que podem fornecer. 

No N. E. E.M. foi desenvolvido um modelo Gaus- 
siano de fontes múltiplas e tem sido realizado trabalho 
conducente a um modelo geral de previsão de poluição 
baseado na resolução das equações da camada limite 
atmosférica e de difusão turbulenta. 


1. INTRODUÇÃO 


A planificação do controle da poluição do ar deve 
ter em conta as relações entre as emissões de substân- 
cias na atmosfera e as concentrações medidas na atmos- 
fera. As repercussões sobre a saúde e os outros efeitos 
nocivos devidos à poluição do ar são proporcionais às 
concentrações de substâncias tóxicas transportadas pela 
atmosfera; os níveis de qualidade do ar ambiente são por- 
tanto geralmente expressos sob a forma de concentrações 
máximas permitidas. Contudo os programas de controle 
só indirectamente podem afectar as concentrações redu- 
zindo as emissões para a atmosfera. 

A avaliação da eficácia dum programa de redução 
das emissões necessita da utilização dum modelo que 
permita transformar os dados relativos às emissões — es- 
pecialmente localizadas —- em dados sobre a qualidade 
do ar em toda a região tratada, 
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ABSTRACT 


This paper is mostly intended to potential users of 
air pollution models and to those more interested in sol- 
ving a particular problem than on the physical processes 
involved. 

lt describes in broad terms the more common use 
of pollution models name!y planification of land use, 
survey on control strategies and alert systems. 

Gaussian models are particularly valuable for sol- 
ving this class of problems as, in spite of their limita- 
tions, they lead to adequate results, are simple to use and 
require little computer time. Reference is made to the 
input data requirements and to the type of information 
the model can provide. 

A multiple source Gaussian model was developed 
at NEEM. Work on a general model based on the solu- 
tion of the atmospheric boundary layer and turbulent dif- 
fusion equations is under way. 


Dum modo geral um modelo de previsão da quall- 
dado do ar pode ser definido como um esquema de cál- 
culo utilizando as informações respeitantes às emissões, 
parâmetros meteorológicos, topográficos, etc., fornecendo 
perfis de concentrações dos vários poluentes atmosféri- 
cos. À maior parte dos modelos operacionais foi conce- 
bida considerando os processos de transporte e de difu- 
são das substâncias na atmosfera e são geralmente cha- 
mados modelos de difusão atmosférica. 

A utilização dos modelos de previsão de poluição do 
ar é geralmente classificada em três categorias: 


1. Planificação para utilização do solo 
Compreendendo as zonas industriais, estradas, 


«cinturas verdes», etc., para obter uma qualidade de ar 
aceitável em toda a região em causa. 


2. Avaliação da estratégia de controle 


Em sentido absoluto (tal estratégia dará uma qua- 
lidade de ar aceitável?) e em sentido relativo (por ex. 
eficácia económica de várias estatégias). 


3. Funcionamento dum sistema de alerta 


Identificando as ocasiões em que são necessárias 
reduções de emissão a curto prazo, 


É necessário ter em conta a precisão que é possível 
obter com a utilização dum modelo de previsão. 

A incerteza de cada tipo de dados (parâmetros de 
emissão, factores meteorológicos e medidas da qualidade 
do ar), associada à incerteza na descrição dos fenómenos 
de transporte atmosférico, leva à adopção de grandes 
margens de erro para a maior parte dos cálculos dos 
modelos. 

Tanto no caso de modelos de previsão a longo prazo 
(planificação, avaliação de estratégias) como a curto 
prazo (sistemas de alerta), a precisão das previsões 
nunca pode estar contida em limites estreitos como por 
exemplo + 10 %. 

A precisão das previsões por modelos é função 
das escalas de tempo e de espaço e da precisão dos da- 
dos fornecidos para a utilização do modelo. Por exemplo, 
as previsões referentes às concentrações médias anuais, 
estabelecidas para uma grande região urbana sobre a qual 
se possuem informações aceitáveis sob o aspecto meteo- 
rológico, topográfico e sobre as emissões, podem apre- 
sentar uma falta de precisão sistemática indo nos mode- 
los correntes até + 25 %. Noutro caso, as previsões rela- 
tivas às concentrações máximas quotidianas no interior 
duma região geográfica mais pequena podem apresentar 
um factor de imprecisão de 2 ou mais. Estes números 
referem-se a modelos de tipo gaussiano. 

Apesar dos elevados erros próprios dos modelos de 
previsão, o emprego destes modelos é a solução mais 
racional para a avaliação dos programas de redução das 
emissões. 


2. DESCRIÇÃO SUMÁRIA DO MODELO GAUSSIANO 
2.1 Modelo Matemático 


Os modelos para a previsão da poluição atmosfé- 
rica mais frequentemente utilizados são baseados numa 
representação empírica das plumas provenientes de fon- 
tes de poluentes. 

Para uma fonte pontual contínua a equação empírica 
da pluma mais frequentemente utilizada em cálculos de 
modelos é o modelo de Pasquill-Gifford. 


c( H) = = ad 
id = 
(1) 
1 z+H 1 z— MH 
exp [——( 2] +exp [——( /] 
2 o 2 o 
z VA 
26 


em que ox coincide com a direcção do vento; oz com a 
direcção vertical e oy transversal. 


04 + € — desvio padrão de distribuição de po- 
luente nas direcções horizontal e verti- 
cal, respectivamente, função da distân- 
cia à fonte x e da classe de estabili- 
dade de Pasquill (m). 

admitindo: 


— velocidade do vento constante u (m/s) 


— fonte pontual de intensidade OQ (mg/s) à altura 
efectiva H (m) 


—. gradiente de concentração nulo ao nível do solo 
dc 


Para as grandes zonas urbanas em que o número 
de fontes individuais é demasiado grande para que cada 
uma seja tratada separadamente é normal o agrupamento 
das fontes de importância relativamente pequena (caldei- 
ras de aquecimento, automóveis, etc.) numa fonte super- 
ficial ou linear, tratando-se individualmente apenas as 
grandes fontes (centrais térmicas, grandes indústrias). 

A extensão do modelo de fontes pontuais de modo 
a incluir também fontes superficiais e/ou lineares é feita 
seguindo a metodologia utilizada na dedução do modelo 
de fontes pontuais, 


2.2 Restrições 
A equação (1) apresenta alguns aspectos essenciais: 


a) À equação representa uma conservação de massa 
e um perfil derivados dum grande número de 
referências «in situ» sobre a dispersão atmosfé- 
rica. 


b) A equação da pluma representa concentrações 
médias; um tempo médio de 1 h pode ser consi- 
derado como representativo, dependendo da al- 
tura da emissão H e da estabilidade das carac- 
terísticas do vento. 


c) A equação é representativa de condições estacio- 
nárias sendo muito menos válida no caso de mu 
danças rápidas das condições meteorológicas. 


d) Supõe-se que o eixo da pluma é horizontal e des- 
preza-se o perfil vertical do vento. 


e) O modelo não toma em conta as perdas de ma- 
terial por lavagem por chuvas, precipitação, 
reacção química, etc. Este factor torna-se mais 
importante à medida que aumenta o tempo de 
permanência na atmosfera e para efeitos a longa 
distância. 


f) O modelo não tem em conta os efeitos de topo- 
grafia local nem da proximidade de vastas su- 
perfícies líquidas. 
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g) A equação não pode ser aplicada a condições 
de calmaria em que u = O. Este é, de resto, o 
principal defeito inerente aos modelos da fonte 
pontual efectiva para aplicação a sistemas de 
alerta; o risco de episódios de forte poluição 
atmosférica é tanto maior quanto mais fracos 
e variáveis são os ventos. 


2.3 Tipos de poluentes utilizados para as previsões 


O emprego de SO; e/ou CO como medida da polui- 
ção não significa que sejam estes os poluentes mais peri- 
gosos do ponto de vista da saúde ou económico mas as 
informações de que se dispõem sobre estas duas substân- 
cias são mais numerosas do que sobre os outros po- 
luentes. 

A maior parte dos modelos de previsão, sejam mo- 
delos a curto ou a longo prazo, têm sido baseados 
nas emissões e concentrações de anidrido sulfuroso. 
O SO;, além de ser um dos poluentes mais abundante- 
mente emitidos, serve como indicador das concentrações 
dos outros poluentes usuais sobretudo nos períodos de 
estagnação atmosférica. 

Quando se quer avaliar a poluição causada pelos 
veículos automóveis é costume utilizar o monóxido de 
carbono que é emitido quase exclusivamente nos gases 
do escape dos motores de combustão interna (entre 80 
e 90 %). 


Quando o estudo diz respeito a áreas essencialmente 
residenciais, comerciais e industriais e/ou centrais tér- 
micas, é preferível escolher o SOs. 

Em relação ao CO, os outros poluentes apresentam 
problemas especiais. O SO, está sujeito a várias trans- 
formações químicas. As partículas atmosféricas provêm 
do muitas fontes naturais e artificiais. 


Quando a outros poluentes emitidos pelos veículos 
automóveis como NO, e hidrocarbonetos, as suas con- 
centrações têm sido medidas em poucos locais. 

Outros poluentes como o chumbo e o arsénio podem 
pôr problemas sérios em algumas regiões mas não têm 
sido utilizados como base para modelos de planificação 
ou sistemas de alerta. 


2.4. Dados necessários ao funcionamento dum modelo 
Gaussiano 


Para calcular a concentração provável num receptor 
devida a uma única fonte, é indispensável ter informa- 
ções sobre a potência e elevação da pluma, as posições 
da fonte e do receptor em relação à direcção média do 
vento e os parâmetros meteorológicos usuais. Para conhe- 
cer o efeito de fontes múltiplas sobre um receptor único 
(caso dum modelo urbano) basta adicionar as contribui- 
ções de cada fonte simples. 

Finalmente, para obter concentrações médias num 
período longo, calculam-se as concentrações possíveis 
resultando dum certo número de condições meteoroló- 
gicas standard e calcula-se uma concentração média 
ponderada. 
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2.4.1 Parâmetros das emisões 


Às técnicas de colheita aos dados das emissões su- 
põem o emprego de: 


a) questionários submetidos directamente aos prin- 
cipais emissores ou mesmo a todas as indús- 
trias implementadas na região controlada; 


b) inspecções e controle das fontes, sobretudo no 
caso dos maiores emissores simples (poucas 
chaminés de grande potência); 


c) recenceamento da população, das habitações par- 
ticulares e das construções comerciais; 


d) recenceamento do tráfego automóvel; 


e) registos industriais que frequentemente catalo- 
gam os processos industriais por tipo, dimensão 
e localização; 


f) dados sobre a produção de energia eléctrica em 

centrais térmicas. 

O tipo de fonte está relacionado com o poluente uti- 
lizado no modelo. No caso do anidrido sulfuroso apresen- 
tam-se normalmente fontes pontuais e superficiais; estas 
substituem o grande número de pequenas fontes ind'vi- 
duais. Se se trata do óxido de carbono, consideram-se 
fontes lineares pois a maior parte do CO provém de mo- 
tores de combustão interna. Se se escolhem outros po- 
luentes como representativos para uma dada área (No, . 
hidrocarbonetos) determinam-se as fontes representativas 
considerando os tipos de fontes existentes na região es- 
tudada. 

No caso de grandes fontes fixas, tratadas individual- 
mente nos modelos de fonte pontual efectiva é desejável 
que se reunam as seguintes informações sobre as emis- 
sões: 

a) Força média da emissão (para cada chaminé ou 

cada processo) 


b) Carácter das emissões (contínuas, interrompidas, 
ocasionais) 


c) Regime típico das emissões no tempo (variabi- 
lidade em função da hora do dia, do dia da se- 
mana, da estação do ano, da temperatura, da 
procura de produção) 


d) Detalhes das emissões (altura e diâmetro da 
chaminé, velocidade e temperatura do efluente, 
outras chaminés e construções próximas) 


e) Reduções possíveis das emissões em caso de 
urgência (combutíveis depurados, paragem pa!- 
cial do processo) 

Este último dado é particularmente útil para os mo- 
delos de alerta; estes dados permitem examinar os me- 
lhoramentos possíveis resultantes das diferentes redu- 
ções de emissão. Podem ser igualmente avaliados os 
custos dos sistemas de alerta e as despesas associadas 
aos diferentes graus de redução das emissões. 

No caso de fontes móveis, um inventário completo 
das emissões destinado ao cálculo por modelos deve 
incluir: 


a) densidade do tráfego, agrupado por tipo de veí- 
culo 


A 


b) taxas de emissões, agrupadas por tipo de vel- 
culo e por ciclo de operação (aceleração, cru- 
zeiro, desaceleração, ralenti) 


c) Ciclos operatórios para cada tipo de veículo, 
para cada ramificação e cruzamento de estradas 
na região estudada. Os dados deverão incluir a 
variabilidade segundo a hora do dia, o dia de 
semana, a estação do ano e as situações espe- 
ciais. 

O inquérito deve ser feito segundo uma forma adap- 
tada ao funcionamento do modelo dependendo da utili- 
zação do modelo (planeamento, alerta, etc.) o detalhe 
nos dados. 


2.4.2. Dados meteorológicos 


Os dados meteorológicos utilizados nos modelos 
gaussianos são de quatro tipos: velocidade e direcção do 
vento, estabilidade da atmosfera e profundidade da ca- 
mada de mistura. 


Em geral os dados meteorológicos são medidos em 
poucos locais, muitas vezes apenas no aeroporto que 
serve a região tratada quase sempre em zona não urbana 
e normalmente plana. O emprego de tais dados meteoro- 
lógicos nos cálculos das concentrações de poluentes em 
zonas urbanas pode introduzir erros importantes cujo 
peso deve ser avaliado para cada região em estudo. 


É desejável que se disponha de dados meteoroló- 
gicos referentes à zona estudada em vários locais, se a 
zona for grande e heterogénea. 


a) Velocidade e direcção do vento 


Na maior parte dos cálculos de modelos, utilizam- 
-so valores únicos para a velocidade e direcção do vento 
na região estudada, muito embora possa haver mais que 
um local de medição. Neste caso, os valores adoptados 
são a média das observações realizadas nos vários locais. 


A velocidade e direcção do vento deveriam ser 
determinadas à altura do eixo da pluma; o mais vulgar 
é serem utilizados valores de superfície visto que são 
raras as medidas aerológicas suficientemente detalhadas. 


Os cálculos da concentração de poluentes por mo- 
delos gaussianos são especialmente sensíveis aos erros 
na direcção do vento pelo que é desejável que se dis- 
ponha de observações da direcção do vento feitas nas 
zonas estudadas. 


b) Estabilidade atmosférica 


Os cálculos dos modelos gaussianos de fontes múl- 
tiplas exigem o conhecimento dos coeficientes de dis- 
persão da pluma 980, 


Estes parâmetros podem ser avaliados a partir de 
medidas da flutuação dos ventos horizontais e verticais 
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ou dos dados de estabilidade atmosférica (método ay 
Pasquill-Gifford) [1]. 


Baseando-se num conjunto extenso de observações 
Pasquill pode relacionar os coeficientes de dispersão com 
condições meteorológicas locais que agrupa em seis 
classes de estabilidade, de extremamente instável (classs 
A) a moderadamente estável (classe F), a partir do va- 
lor da velocidada do vento à superfície e da insolação 
diurna ou cobertura de nuvens nocturna. 


Este método não é tão fiável como as determinações 
a partir da sondagem de temperatura; em contrapartida, 
é de fácil e rápida aplicação, sendo estes dados msteo- 
rológicos disponíveis em muitas zonas. 


A medida que se vão instalando equipamentos de 
medida dos ventos para dar apoio a programas de fisca- 
lização de poluição atmosférica, é desejável que se mon- 
tem instalações e serviços capazes de medir e registar 
continuamente a flutuação dos ventos. 


c) Profundidade da camada de mistura. Altura da 
base da camada de inversão 


Os dados sobre a profundidade da camada de mis- 
tura são obtidos a partir de sondagens verticais de tem- 
peratura e da temperatura superficial. 


As sondagens verticais de temperatura para aplica- 
ção em fenómenos de dispersão de poluentes devem ser 
feitas com pequenos incrementos de altitude na zona da 
baixa camada limite atmosférica (primeiras centenas de 
metros) de modo a que se obtenha um perfil com sufi- 
ciente resolução. 


Quando não existem dados de sondagem vertical 
numa dada região, podem utilizar-se dados provenientes 
de estações de observação distantes para avaliar a pro- 
fundidade da camada de mistura (altura base da ca- 
mada de inversão). O perfil de temperaturas da atmos- 
fera a partir de altitudes de cerca de 200 m passa a ser 
geralmente controlado por fenómenos meteorológicos à 
escala sinoptica. É portanto possível fazer interpolações 
aceitáveis de perfis de temperatura vertical para distã- 
cias horizontais podendo atingir 300 km. No entanto a 
utilidade destes dados extrapolados é limitada por duas 
influências: 


a) a estrutura da atmosfera inferior pode ser modifi- 
cada por efeitos de grandes aglomerações urba- 
nas, 


b) quando as condições atmosféricas estão a mudar 
as sondagens verticais feitas a 300 km da região 
estudada podem não ser representativas desta 
região. 


2.4.3. Topografia 


As características topográficas influenciam muito 
os fenómenos de transporte no interior da camada limite 
planetária. Os modelos gaussianos de fontes múlt'plas 
não têm em conta os efeitos da topografia nem qualita- 
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QUADRO | 


ESQUEMA DOS PLANOS DE CÁLCULO PARA MODELOS DE «FONTE PONTUAL EFECTIVA» 


|. PARA CADA PAR 
FONTE-RECEPTOR 


PARÂMETROS DA FONTE PARÂMETROS DO RECEPTOR PARÂMEROS METEOROLÓGICOS 
1. Taxa de emissão. . 1 Localização. | 1. Direcção e velocidade do 
: : vento. 
2. Altura efectiva da fonte. a. Distância à fonte. 
b. Localização em relação 2. Classe de estabilidade. 


3. Localização efectiva da 


3. Profundidade da camada 


à direcção do vento 
fonte. 


2. (Durarão da exposição). de mistura. 


, 


EQUAÇÃO DA PLUMA NO | 
ESTADO ESTACIONÁRIO 


Concentração possível no recep- 


tor devida à fonte 


it, PARA VÁRIAS FONTES 


AFECTANDO O RECEPTOR Adicionar a contribuição de cada fonte individual. 


Calcular as concentrações correspondentes a cada elemento num grupo standardizado de direcções de vento, velocidades 


HPARA: UMA CONCENTRAÇÃO de vento e medidas de estabilidade atmosférica. 


MÉDIA DURANTE UM PERÍODO 


LONGO, NO RECEPTOR, Em seguida calcular a concentração média ponderada pela frequência normalizada de ocorrência para cada uma das 
DEVIDA A VÁRIAS FONTES combinações dos elementos meteorológicos. 
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tiva nem quantitativamente, mas estes efeitos podem 
ser tratados por meios «ad hoc» quando se dispõe de 
medidas «in situ» para verificar e melhorar os cálculos 
dos modelos. 


2.4.4. Medidas da qualidade de ar 


Os dados relativos às medidas da qualidade do ar 
não sendo necessários como dados de entrada (condi: 
ções iniciais) para o funcionamento do modelo de fontes 
múltiplas, são indispensáveis ao melhoramento e verifica- 
ção dos esquemas de cálculo do modelo. A variabilidade 
dos fenómenos locais, afectando os processos de trans- 
porte e de dispersão atmosféricos é tão ampla qua se 
torna necessário testar, para cada região, as possibilida- 
des dos modelos. 


As medidas das concentrações de poluentes apre- 
sentam bastantes dificuldades. Mesmo quando os apare- 
lhos funcionam com uma precisão óptima, 2 apare'hos 
iguais colocados lado a lado numa estação de observa- 
ção fornecem muitas vezes valores diferentes, À correla: 
ção entre os resultados dos modelos e as medidis da 
qualidada do ar não podem ser melhores que o grau de 
reprodutibilidade das próprias medidas. 


3. EMPREGO DO MODELO DE POLUIÇÃO DO AR NO 
ESTABELECIMENTO DE SISTEMAS DE ALESTA, 
A AVALIAÇÃO DA EFICÁCIA/CUSTO DAS ESTRA- 
TÉGIAS DE CONTROLE 


Dadas as repercuções na saúde pública devidas à 
exposição das populações e também plantas e mater.ais 
a concentrações elevadas de poluentes, é indispensável 
que os organismos de controle possam prever esses pe- 
ríodos a fim de reduzir as emissões na expectativa de 
condições desfavoráveis. 


Um dos aspectos em que os modelos de previsão 
da poluição do ar são particularmente úteis é na prepara- 
ção de previsões a curto prazo e estabelecimento da sis- 
temas do alerta em regiões urbanas. Essas previsãas re- 
ferem-se a situações esperadas dentro de um ou dois 
dias e são geralmente consideradas como «previsões 
para 36 horas». 


3.1. Precisão das previsões 


Para explorar ao máximo as possibilidades dum mo- 
delo do previsão, reduzindo ao minimo os riscos de de- 
cisões erradas, a avaliação da precisão deve fazer parte 
das informações que servem de base às decisões relati- 
vas aos controles de urgência, Uma análise regular da 
precisão resultante de previsões correntes é o melhor 
método de avaliação da precisão dos métodos de pre- 
visão. 

As correcções podem ser propostas a introduzidas 
pelo pessoal que opera o modelo o qual fica a connecer 
as possibilidades e limites do mesmo e pode formular 
os melhoramentos empíricos às técnicas de previsão. 
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QUADRO 


GAMAS DE ERROS TÍPICOS PARA OS MODELOS DE 
POLUIÇÃO ATMOSFÉRICA DO TIPO 
«FONTE PONTUAL EFECTIVA» 


E 


Estações Médias 
individuais urbanas 
Cálculos a 
curto prazo 1,0 0,5 
(de 1 hora a 1 dia) 
Cálculos a 
longo termo 0,35 0,2 


(de 1 mês a 1 ano) 


A gama do erro é definida como a média RMS de 
conc. observada — conc. calculada 


conc. observada 


Como já foi referido e se conclui do Quadro li, a 
precisão dum modelo aumenta como periodo de tempo 
a quo ele é aplicado e com o tamanho da região. 


3.2 Informações a obter dos cálculos dos modelos 


Geralmente, os cálculos de modelos dão as redes 
possiveis do concentração na região controlada deste o 
instanto presente até vários instantes futuros. 

O modelo de previsão pode fornecer os seguintes 
tipos de informação específicas: 


1) mapas indicando os sectores da região contro- 
lada em que as normas aplicáveis à qualidade 
do ar serão provavelmente ultrapassadas; 


2) valores previstos para as concentrações máximas 
e suas localizações; 


3) concentrações previstas nos centros povoados 
da região controlada; 


4) duração prevista do período de forte concentra- 
ção; 
b) redes de concentração resultantes dos vários 
controles de emissão possiveis. 
Na maior parte das vezes as informações necassá- 
rias são fornecidas mais eficazmente na forma de mapas 
preparados em computador. 


3.3 Período de tempo necessário para desencadear um 
alerta 


Os períodos de tempo necessários para desencaasar 
um alerta são determinados pelo «tempo de reacção» 
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para estabelecer as medidas de controle e são limitados 
pela duração das previsões meteorológicas a utilizar 
como dados de base para o cálculo dos modelos. 

O tempo de reacção para alertas a curto prazo pode 
ser de algumas horas ou durar um ou dois dias, segundo 
os tipos de emissões a controlar, a existência dum plano 
de urgência pré-estabelscido e as comunicações possí- 
veis. 

O tempo total de reacção é a soma de três parcelas: 


a) Duração das comunicações 


Se so trata de alertar apenas algumas grandes fon- 
tes (centrais térmicas, grandes indústrias), o tempo é 
muito curto. À medida que aumenta o número de fontes, 
aumenta o tempo de comunicação; no caso de toda uma 
população (redução do tráfego automóvel) a transmissão 
pode durar meio dia. 


b) Tempo de execução 


O tempo necessário para a execução de ordens de 
controla das emissões pode variar entre uma e várias 
horas, de acordo com as operações necessárias (modifi- 
cação de combustível, interrupção de processo, paragem 
de tráfego), etc. 


c) Tempo necessário para que a redução das emis- 
cões afecte a qualidade do ar. 


Como os períodos de forte poluição provocada por 
emissões a baixa altitude correspondem a períodos de 
vento fraco, podo haver várias horas de diferença entre 
as reduções das emissões e uma melhoria notável de 
qualidade do ar local. 

Determina-se um «tempo de ventilação» para uma 
determinada zona dividindo um comprimento caracteris- 
tico da zona pela velocidade média do vento. 


3.4. Emprego de modelos de sistema de alerta para a 
avaliação da eficácia e custo do controle 


A muio: parte das medidas urgentes de controle a 
curto prazo que os organismos podem impôr (por ex., in- 
terrupção da produção industrial e dos transpotes, em- 
prego de combustíveis depurados) introduzem encargos 
cuplemontares à economia, Naturalmente as decisões de 
impôr tais controles devem apoiar-se sobre uma base 
adequada de dados e conhecimentos experimentais. 

Dispondo um organismo dum modelo de controle 
para a previsão da poluição a curto prazo, modificado 
e verificado para a região, o organismo pode utilizar esse 
mesmo modelo para avaliar a eficácia das diferentes me- 
didas de controle de urgência propostas para a região. 

Também estes modelos podem ser utilizados para 
estimar o custo das medidas de controle de urgência a 
adoptar desde que se conheçam as várias opções de con- 
trole e se tenha estimado o preço de todos os elementos 
deste plano de controle: preço de combustíveis depura- 
dos para as centrais ,refinarias e instalações de aqueci- 
mento de grandes edifícios comerciais, perda de produ- 
ção resultante de abrandamento ou interrupção duma ope- 
ração industrial, o valor do tempo perdido se o pessoal 
de empresas comerciais e industriais largar o trabalho 
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antes da hora de saída e o velor da perda do trabalho nos 
centros comerciais cujas actividades são entravadas pelas 
restrições obrigatórias nos transportes. O modelo pode- 
ria ser utilizado para avaliar o custo dos diferentes planos 
e apesentar estes dados aos responsáveis dos organismos 
de controle que, bem informados, poderão tomar a de- 
cisão mais adequada considerando ainda factores que 
não foram incluídos quantitativamente na análise feita 
pelo computador. 


4. NECESSIDADE PARA O DESENVOLVIMENTO E EM- 
PREGO DE MODELOS 


4.1. Normas de qualidade de ar 


É necessário que sejam adoptadas normas de quali- 
dade de ar ou valores indicativos recomendados a ser 
obedecidos pela atmosfera da região a modelar. 


4.2. Implementação do modelo gaussiano em computador 


Foi desenvolvido no N. E. E. M. «Um modelo gaus- 
siano para a dispersão de poluentes na atmosfera no 
caso de fontes múltiplas pontuais e inversão de tempera- 
tura» [2]. No modelo estabelece-se uma rede horizonial 
uniforme sobre um mapa da região a estudar. As malhas 
são quadradas e nos modelos correntes a dimensão do 
lado da malha tem variado entre 0,5 e 10 km, segundo 
a aplicação do modelo e a precisão dos dados fornecidos 
ao modelo. Malhas de 1 e 2 km têm sido usadas com 
frequência. 

O modelo [2] está implementado num calculador 
Hewllet Packard 30 acoplado ao qual existe um traçado 
de gráficos HP 9862 que permite que as saídas se façam 
sob a forma de mapas de linhas de isoconcentração. 
Esta forma é geralmente utilizada pois permite uma ve- 
rificação rápida de quais as zonas atingidas por altas 
concentrações. 

O calculador HP30 é programável em BASIC pelo 
que a implementação do modelo em FORTRAN — quase 
universalmente usada em cálculo científico — é simples 
o rápida. O calculador tem 5856 palavras de memória e, 
na fase actual, o modelo admite 10 fontes pontuais e 
1296 nós de rede O número de fontes pode ser aumen- 
tado diminuindo o número de malhas da rede e vice-versa 
para o mesmo tempo de cálculo. 

Sendo necessário utilizar um número de fontes su- 
perior, é possível fazê-lo sobrecarregando o tempo de 
execução no computador. 

Como foi acentuado, este modelo pertence a um tipo 
de modelos bastante grosseiros com um número apreciá- 
vel do restricções. No entanto têm sido e continuam a 
ser extensivamente utilizados nos U.S. A. e U.K. como 
base a sistemas de alerta, para planificação de regiões, 
e avaliação de estratégias de controle. 

No N. E. E. M. têm sido desenvolvidos métodos nu- 
méricos de diferenças finitas que permitem a resolução 
de sistemas de equações às derivadas parciais de tipo pa- 
rabólico [3] a [6] assim como modelos de turbulência 


[7]. 
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Estes meios permitem-nos resolver as equações fun- 
damentais do escoamento da atmosfera na zona em que 
so dispersam os poluentes (camada limite), assim como 
a difusão turbulenta desses memos poluentes. 

A vasta literatura publicada sobre o assunto indica 
claramente uma evolução, dos modelos simplificados do 
tipo gaussiano ou outros, para os métodos baseados na 
resolução das equações fundamentais. No entanto estes 
métodos são morosos, exigem computadores de grande 
capacidade, dados muito detalhados e precisos quanto às 
taxas de emissão e condições meteorológicas. 

Os modelos deste tipo encontram-se ainda em fase 
de investigação e desenvolvimento. 


4. 3. Dados necessários 


Os dados necessários ao funcionamento do modelo 
gaussiano podem classificar-se em: 


a) Parâmetros das fontes —- sua localização e ca- 
racterísticas, nomeadamente: taxa de emissão, 
altura, diâmetro interior da boca, velocidade de 
saída e temperatura dos gases à saída. 


b) Parâmetros meteorológicos — velocidade e di- 
recção do vento, classe de estabilidade segundo 
Pasquill, altura da camada de mistura ou altura 
da base da camada de inversão e temperatura 
da atmosfera. 


c) Medidas da qualidade de ar — não sendo neces- 
sárias ao funcionamento do modelo, são-nos 
para teste, validação e ajustamento do modelo 
à região a que é aplicado. 

Qualquer que seja o tipo de modelo de previsão que 
se utilize, os dados referidos são sempre necessários po- 
dendo variar um pouco a forma sob que são apresean- 
tados. 
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A recolha dos parâmetros das fontes é uma tarefa 
morosa e pesada; muitas vezes é nesta zona em que é 
necessário investir mais esforço pelo que é aconse- 
lhável iniciar a recolha duma forma sistemática e adap- 
tada ao funcionamento do modelo com a antecedência 
possível. 
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" Electrochemical basis of metallic 


corrosion in melts” 


ABSTRACT 


After pointing out the electrochemical nature of cor- 
roston reactions in molten salts, an account of the essen- 
tial parameters concerning the electrochemical corrosion 
phenomena is presented. Emphasis is given to the elec- 
trodics of corrosion; based on the «Erdey-Gruz-Volmer» 
theory and on the «acid-base concept» of Lux and Flood 
a relationship between the corrosion potential and the 


function pO “is derived. Diagrams in which the equili- 
brium potential E in the metal-fused salt system is plot- 
ted against the activity of oxide ion in the melt are de- 
veloped, using the Pourbaix method of data presentation; 


as an example, an E-pO*  diagram is constructed for nic- 
kel at 900º C, in contact with molten sodium sulphate, 
and its limited nature as an aid to corrosion control is 
summarized. The need for study of metal corrosion in 
terms of electrochemical kinetics is indicated. Further, 
it is suggested to apply the newer electrochemical techni- 
ques to the study of corrosion kinetics in high-tempera- 
ture media. 


ELECTROCHEMICAL NATURE OF ATTACK 


Theoretical considerations of the corrosion pheno- 
mena have often been considered inaccessible, mainly 
due to the great complexity of the corrosion processes. 
It is generally agreed today that several factors with dif- 
ferent physicochemical characteristics can play an impor- 
tant role in the mechanisms of corrosion. However, des- 
pito these difficulties, corrosion problems have been su- 
ccessfully interpreted in recent years, 

Particularly over the last thirty years, corrosion stu- 
dies have demonstrated that the primary process bet- 
ween the corroding metal and corrodent is electrochemi- 
cal in nature; attack is basically a chemical reaction ac- 
companied by the transfer or displacement of electrons 
(Evans and Hoar 1932). Electrochemical theory today 
rests secure in its application to corrosion phenomena in 
aqueous media (Evans 1960). 

If molten salts are regarded as infinitely concentra- 
ted aqueous solutions it is a natural step to try to formu- 
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RESUMO 


Uma vez apontada a natureza electroquimica da 
corrosão metálica em meio fundido, descrevem-se os pa- 
râmetros teóricos essenciais da termodinâmica e cinética 
electroquímica da corrosão. Dá-se relevância à electro- 
cinética da corrosão; com base na teoria de «Erdey- 
-Gruz-Volmer» e nos conceitos de ácido e base segundo 
Lux-Flood, estabelece-se uma relação entre o potencial 
de corrosão e o cologaritmo da actividade oxigeniónica 
do electrólito fundido. Introduzem-se os diagramas de 


equilíbrio E-pO"”” , cuja construção se faz segundo o mé- 
todo de Pourbaix, de representação de dados termodiná- 
micos; exemplifica-se com a construção de um diagrama 


E-pO* para o sistema Ni-Na;SO, a 900º C. Finalmente, 
sugere-se que se encare a corrosão metálica nas misturas 
fundidas, em termos da cinética electroquímica, recor- 
rendo-se às técnicas transientes, usadas em electroaná- 
lise, para o estudo da cinética e mecanismos dos pro- 
cessos corrosivos. 


late a description of corrosion in them on the basis of 
concepts like those which have been developed for aque- 
ous solutions. The main difficulty here is that the pro- 
blem is not completely similar to an aqueous system, 
and it is, therefore necessary to assess other electroche- 
mical parameters. Besides, whilst in aqueous solutions 
metais are virtually insoluble in the electronated state, 
in fused salts they can be appreciably soluble, corrosion 
occurring, in such cases, without the de-electronation of 
tho metal. In spite of this the electrochemical approach 
for molten salt systems was put forward about forteen 
years ago and is finding general acceptance. Experiments 
carried out by Baudo et al. (1970), Burrows and Hills 
(1966, 1970), Cutler (1971), Janz and Cont (1964), 
Littlewood and Argent (1962), Sequeira (1974), et al. 
have contributed considerably in this respect. 

There have been some objections concerning elec- 
trochemical interpretations of corrosion based in the 
thermodynamics and kinetics of electrochemistry. In fact, 
these theories are restricted to cases well defined and 


checked by experience. They cannot, however, include 
all the simulteneous corrosion processes. Following these 
essential considerations, some basic concepts of electro- 
chemistry are discussed; modern electrochemical views 
of corrosion in molten salts are also included and a syn- 
thesis of all the miscellaneous aspects considered have 
been put together. 


THE SINGLE POTENTIAL OF AN ELECTRODE 


In the most general sense of the word, an electrode 
is a system consisting of two phases in contact with 
each other, which can be the seat of an electrode reac- 
tion, |. e. a reaction in which certain constituents of the 
two phases participate and by which a transfer of charge 
takes place from the bulk of one phase to the bulk of the 
other. In a mora restricted sense an electrode is defined 
as a system metal /electrolyte or metal /solution of one or 
several electrolytes. 

it has been shown over the last 70 years that at the 
interfaca of two phases one finds an electrical double 
layer with a characteristic potential (Adam 1938, and 
Butler 1940). The origin as well as the characterisation 
of the electrochemical double layer on a metal immersed 
in an aqueous electrolyte have been extensively studied. 
Similar studies in melts have been carried out only recen- 
tly. This paper is not intended to deal with this aspect in 
detail but the subject is mentioned because of its impor- 
tance. Very good reviews on the subject have been pu- 
blished by Devanathan and Tilak (1965), Graves et al. 
(1966), and Ukshe et al. (1964). 

Returning to the subject of this section, let the inner 
electric potential of the metal be Pm and that at a re- 
mote point in the molten electrolyte be & . Moreover, let 
* + represent the Galvani electric potential across the 
electrochemical double layer, i.e. 4 Pb = 4, — P,. When 
a dynamic equilibrium across the double layer is set up, 
the electrochemical free energy of one mole of metal 
in the surface will be equal to the electrochemical free 
energy of one mole of ions on the melt side of the dou- 
ble layer: — 


um E bs (1) 
or Bo TOP d” nF&. 
and bm —u,=—nFAS (2) 


um — Hg is the chemical energy balanced by the elec- 


4 


trical energy liberated as the ions M"* traverse the 
doublo layer; F is Faraday's constant and n is the num- 


ber of free electrons exchanged at the Mº C4M inter- 
face. It is impossible to measure à & directly but, if the 
metal/melt electrode system is coupled «aback-to-back» 
with a second arbitrarily chosen electrode system, it is 
possible to obtain a relative potential difference. If a zero 
potential and a convenient reference scale of potentials 
is defined, the combined system acts as a cell between 
whose electrodes there exists an e.m.f. given by 


E, 1 =AG-—- AC, (3) 
(assumed zero) 


E, is called the single potential of the M"*/M elec- 
trode on the chosen standard reference scale. Conside- 
ring now a more concise electrode reaction 


Ss; uM +ne” =0 (4) 


where M, represents the constituents of the electrode 
taking part in the reaction and e - represents the electron, 
the corresponding conditions of dynamic equilibrium can 
ba expressed in the general form 

DS. vu; —n E E n4 = O (5) 


i e E 


Relating the chemical potentials, u; , to the standard 
chemical potential, u'; , equation (5) can be put into the 
form 

RTX v logiM,); = —Z »; ptnF Es (6) 
Changing from natural to Briggsian logarithms, the last 
relation can be rewritten in the form 


2.3026 RT 


ad =; vilog iM,;) 


(7) 


in which E Gs has the following value in terms of the 
standard chemical potentials of the reactants 


0.4 n F (8) 


Equation (7) is the Nernst equation. Ejs is the single 
potential of the electrode measured when the activity of 
every constituent of the system is equal to unity, and is 
known as the normal or standard potential, Once the nor- 
mal electrode potential on the standard reference scale 
is known, the Nernst equation permits a calculation of 
what the equilibrium electrode potential will be when the 
system has any (M,) + 1. Extensive analytical treat- 
ments of this subject are given by Conway (1965), De 
Groot (1951), Grahame (1947), Guggenheim (1929), 
Parsons (1954), Prigogine (1947), et al. Experimental 
proof of the validity of the Nernst equation in molten 
media is given by Blander et al. (1959), Duke and Gar- 
finkel (1961), Gordon (1899), Ranford and Flengas 
(1965) and many other workers. 


THE TAFEL RELATIONSHIP 


Considering now the dynamic side of the phenome- 
non, particular importance will be attached to electro- 
chemical kinetic studies. A classic relationship precisely 
relating a characteristic parameter of the electrode, its 
electrochemical potential, and the current density tra- 
versing it, is that found by Tafel in 1905. Thereby this 
semi-empirical law enables us to establish a connection 
between the «thermodynamic electrode potential» and 
the current density which defines the kinetic character 
of the phenomenon. 

The Tafel equation may be written 

E=a+blogi (9) 
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E a a a" 


where E is the relative electric potential of the electrode 
ar a current density | (also called «reaction electrode po- 
tential»), a is a constant characteristic of the electrode, 
and b is the Tafel slope, which is one of the parameters 
indicating the mechanism of the electrode reaction. 

Thus, a definition of overpotential for the reaction 
can be introduced by the expression: 


7=E-E, (10) 


where E, is tho reversible or equilibrium electric poten- 
tial possessed by the electrode at electrochemical equi- 
librium (i.e. at zero imposed current) or, more generally, 
the equilibrium potential of the electrochemical reaction 
of the type (4). Equation (10) defines the overpotential 
and considering (9) it is seen that Tafel's equation is it- 
self closely related to the overpotential (y). These pro- 
blems of electrochemical kinetics are discussed in more 
detail later. 

The brief discussion given above is significant 
enough for the present purpose of showing that the Tafel 
relationship leads to the experimental polarization curves. 
Figures 1 and 2 represent the anodic and cathodic poten- 
tial-current curves for two typical electrochemical reac- 
tions, with either an arithmetic scale for the reaction- 
-curren: (Figs. 1-a and 2-a) or with a logarithmic one 
(Figs. 1-b and 2-b). The exchange current density (i, ). 
which is defined as the exchange rate per unit area of 
the potential-determining electrode process at equili- 
brium, is represented in Figs. 1-b and 2-b. The aspect 
of the polarization curves is very important, because it 
gives an account of the degree of irreversibility of the 
reactions. In fact. in the case of a reversible reaction 
(Fig. 1), the application of very small polarizing poten- 
tials is sufficient for producing significant densities of 
current. In the case of irreversibility (fig. 2), on the con- 
trary, tho existence of high overpotentials is a measure 
of the extent of the irreversibility. The reversible or fast 
electrode processes exhibit high exchange current densi- 
ties; on the other hand, the irreversible or slow ones 
exhibit low exchange current densities. Thus, if an elec- 
trodo is more likely to behave reversibly the higher is 
its intrinsic exchange current density. 


Tha Tafel curves, giving E as a function of log 1, 
also provide interesting results about the kinetics of the 
phenomena. In fact, the plot of E against logi is not usu- 
ally a staight line, i.e. the kinetics of the reaction does 
not obey the Tafel Law. The experimental curves obtained 
can be considered as built up of rectilinear portions, or 
on the contrary as showing systematic departures from 
the linear Tafel law for certain current densities reached. 
It is well established in electrochemistry that the Tafel 
law does not apply precisely for small and large current 
densities. One reason for this is that the Tafel equation 
is derived purely by considering the overvoltage due to 
polarization caused by activation energy requirements, 
which assumes that a cation moves away from the metal- 
-fused electrolyte interface as soon as it is dissolved. This 
is generally not the case. The removal of the anodic pro- 
ducts does not increase in the same proportion as the 
current density and the concentration of these products 
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will increase and cause a back e.m.f. This efect is parti- 
cularly relevant in the case of fast electrode processes 
as observation of the Figures 1 and 2 shows: — The dif- 
ficulty in obtaining a pure Tafel slope, i.e. unaffected by 
double layer and mass transfer effects, increases with 
tha kinetics of the electrode processes (i.e., with high 
io | due to increasing overlapping of the diffusion zones 
over the Tafel one. Polarization measurements also in- 
clude ohmic overpotential which arises from an | R drop 
through a portion of the electrolyte between the test elec- 
trode and the reference one, but this term is not signi- 
ficant in molten media as compared with an aqueous 
system. The validity of Tafel's equation has been verified 
by Laitinen and Gaur (1957), Randles and White (1955), 
et al, for many fused systems, 

The facts mentioned above are very important es- 
pecially due to their repercussion on the interpretation of 
electrochemical phenomena. So it is seen, by this simple 
enumeration, that it is possible to connect thermodynamic 
considerations to kinetic ones. 


GENERAL CONSIDERATIONS OF ELECTRODICS OF 
CORROSION 


in section above, when some concepts of electrode 
kinetics where discussed, a quantity of crucial importance 
to electrodics was defined: the overpotential (y). In elec- 
trochemical kinetics the only reaction directly affected 
by the potential is the charge-transfer reaction, i.e. the 
reaction in which charge carriers are transferred across 
the electochemical double layer at a phase boundary. The 
rate of the charge-transfer reaction determines the char- 
ge-transfer, or activation overpotential. This kind of over- 
potential has been treated in detail by Audubert (1942), 
Bockris (1954), Butler (1924), Conway (1965), Erdey- 
-Gruz and Volmer (1930), Gerischer (1960), et al. 

Prior to discussion it is noted that only the basic 
electrodic equations will be considered. Considering the 
conversion of M+* ions into metallic M, the number of 
moles of positive ions reacting per second by crossing 
unit area of the melt/metal interface is proportional to 
the electronation-current density i — and as a first appro- 
ximation (Erdey-Gruz and Volmer 1930): 


| o=—K' (Mt) exp (—-46* /RT) 
=—K' (Mt) exp (— 4G*t / RT) exp — 
- 7=-a FE/RT) 

=—K. (Mt) exp [—(1-a)FE/RT] (11) 


A G* is the activation energy required for the electrona- 
tion | process at the potential E, K' andK are constants 
(M+) is the activity of cations M+ on the melt side of 
the interface, AG* is the chemical activation free 
energy barrier when there is zero electric field acting on 
the ion, and a is the thansfer coefficient for the de-elec- 
tronation reaction introduced by Erdey-Gruz and Volmer 
(1930) which has values ranging from zero to unity 
(0 <a< 1) (Vetter 1952). K is dependant on the re- 
ference potential chosen since K = K exp (á 2) / 


Electronation|De-electronation 


Fig. 1 —- Polarization curves of a fast electrochemical reaction 


Fig. 2 — Polarization curves of a slow electrochemical reaction 
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/ RT). Moreover, positive ions can move from the melt 
to the electrode and also in the opposite direction. Thus, 
there is an electronation reaction 


Mt+ + e -—>M 
and also a de-electronation reaction 
M -—> M+t + em 


In this case, the total current density passing through the 
double layer is 


i=i, +i. 


E Nu (IM) exp (FE /RT) —-K. MT ap [—(l — 
- ar RAR TI (12) 


(M) is the activity of rhe metallic atoms immediately at 
tha electrode surface (unity). If the metal-melt interface 
reaches its equilibrium state there is no net current 
(i= O) but there is an exchange current. The electrona- 
tion and de-electronation reactions continue to occur-but 
at the same rate. The current densities corresponding to 
these individual reactions become the exchange current 
density a (is = N =|i 1). Hence, 


io 1, =K, IM) exp («FE /RT) 
— (Mt) exp [— (1 —a) FEG/RT] (13) 


where E , is the potential difference across the interface 
a* equilibrium. On dividing equation (12) by the corres- 
ponding expressions in equation (13) and taking into 
consideration the definition for overpotential y = E — 
— E, the following important relationship is obtained 
(Erdey-Gruz and Volmer 1930): 


otexp laFmn/RT)-expl—(1—a) Fn/RT]) 
(14) 


For higher charge-transfer overpotentials In | >D> RT /F 
the first or second term on the r.h.s. of equation (14) can 
be neglected depending on the sign of the current so that 
if this equation (14) is put into a logarithmic form a li- 
near relationship results between 7 and log iil similar 
to the Tafel equation: 


2.3026 RT 2.3026 RT 
dg S = E log 1 + , log Ly 
or (15) 
2.3026 RT 2.3026 RT 
E e DO qe, mr MD 
T egr tE g=gF 


Here the subscripts «+» and «-» to the overpotential 
denote anodic and cathodic polarisation respectively. 


APPLICATION TO MOLTEN SULPHATES 


Let us now consider the following reaction of a me- 
tal M in a molten sulphate (Burrows and Hills 1966): 
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1 
M + 7 SO; - Mt + SU o 


(16) 


This overall reaction can be considered as the sum of the 
two half reactions: 


M= Mt + em (17) 


1 
dae= >= 80:* — O 
: 4 E | (18) 


1 

-— "SO, 

2 Pra 
For these two processes to take place simultaneously, it 
is necessary and sufficient that the potential difference 
across the interface be more positive than the equilibrium 
potential for the reaction (17) 


2.3 RT 


+ —>— log (Mt) 
E g 


and more negative than the equilibrium potential of the 
electronation reaction (18) involving electron acceptors 
(such as SO.) contained in the melt 


pay BAN 
os" “08 "sp leg Pso;/Pso, — 
2.3 RT , 
log 10º" 
SF g 


The behaviour for the potential-current curve E (1) which 
can be envisaged for the process (16) is illustrated in 
Fig. 3 by the continuous thick line. | (E) is composed 
additively of the metal-dissolution current |;; (E) and the 
electronation current |; (E), and the two are independent 
of each other at the same potential value. This principle 
of additive combination of all partial processes at an 
electrode surface to obtain the total potential-current 
curve was first formulated by Wagner and Traud (1938). 
The potential E at zero current is neither Fo, 17 Nor E,, 8 
but a so-called mixed potential (E ) for which the total 
potential-current curve gives | (Em ) = O, The total me- 
tal-dissolution current |,; and ciacitodiáiion current l|,; are 
then equal in magnitude but opposite in sign, i.e., 


The rate of corrosion of the metal is obviously given 
directiy by the rate of metal dissolution; hence, the free 
corrosion current |,. (i.e., the rate at which a metal des- 
troys itself) is given by 


Hé = gym — kg 


Considering the electrode kinetics, we can thus write: 


— = — + 
e ” A, E A, Ur Eat 


- A. K, (Mex (a FE./RT)— 
-AK. (Mtexpi-(1-a )JFEÇ/RT] (19) 
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dé: a a (2) 
Egas f--- 
a' 

: A É 

Mº Tc SA. 

a 7 b 

| | ny |piesecuseses 7 

o 

| | ro 
Í log 1 

1 LEO s Ts is WS + bo 
7 Aba tg 7 0 + eta] M Mt + e- (b) 
á | 10 E 4 o at , 
> 50,+.€ mito 550, E! (q!) M o M (b'] 


Fig. 3Corrosion of a metal M in a molten sulphate, 


with 30, -evolution. Anodic and cathodic polari- 

sation curves. Corrosion-potential and corrosion- 

current. 

(1) Representation with an arithmetic scale for 
the reaction- current. 

(2) Representation with a logarithmic scale for 


the reaction- current. 
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= 2 ga 
= 1, =—A, Kas Po, 197 


exp (a F Ec! RT) + 


a 


| 
o 


where As; is the sink area and As is the electron-source 
area (at the interface metal-electrolyte). Assuming that 
overpotentials are sufficiently large, equations (19) and 
(20) become 


— A. K+a7 (M) exp (a, F E. / R T) 


l 


2 


+ is 15 


— q JFEÇ/RTI] 


The free corrosion potential, E E .e. the uniform poten- 
tial difference over the surface of the corroding metal 
can be evaluated from the expressions above; it follows 
that: 


23RT As Ko 
REU -RIÍIP O 4 ra 
La,,+ 1 2) ] 17 K pa UM) 
23RT RA 
fis E Mim 
Clatti-ad]F * sos (21) 


It is clear from this that the free corrosion potential and 
the mixed potential are identical. 

In accordance with the theory of Lux (1939) and 
Flood et al. (1952) who assign a definite acid-base equi- 
librium to all oxyanionic melts, the acid-base equilibria 
for sulphates can be expressed by the equation 


sos” = 80, * 0%” (22) 


and by the equilibrium constant 


K.= Pso, (0%) (802) —! (23) 


Do! 
- 


Assuming [SO | =1, it follows that 


| 
"a 1 1 
lo P o = — o K + e O 2 — 
nda = 4 2 ? (24) 
where pO— = — log [0*-) is a nonprotonic function 


which measures the acidity of the melt. The expression 
(24) for the active species SO, can be introduced into 
Eq. (21), giving 


2 
2.3 RT Ass Kg K 
Saque ad = pesto — 8 18 
Ec=la +19 ca )]F Au Km (M) 
2.3RT 
TETE TT cai 
2a + (1a) ]F (25) 


Assuming that tha ratio As/As; as well as a and ay are 
constant during the overall corrosion reaction (which is 
not free from error), equation (25) becomes 


E. ="4:+ b pO- 


a and b being two constants. Thus, we attain a linear 
relationship between the corrosion potential and the fun- 
ction pO*-, which should be valid for all oxyanionic 
melts. 

What has been presented above is a very elementary 
account of electrodics of corrosion under quite ideal con- 
ditions. The details of the complex corrosion phenomena 
in real systems is out of the scope of this review. 


POTENTIAL — pO'— DIAGRAMS 


The major importance of the function pO*— is made 
clear from theoretical treatments such those published 
by Edeleanu and Littlewood (1960), and Littlewood 
(1962), for example. The basis of these studies rests on 
the usefulness of the diagrams resembling those of Pour- 
baix an his collaborators (1966) for corrosion in aqueous 
solutions, with the modification that the basicity of the 
system is written in terms of pO*— instead of pH. The 
diagrams E—pO'- give information only on the thermo- 
dynamic feasibility of electrode reactions. It follows that 
it is not possible to draw conclusions from the diagrams 
on whether slow kinetic steps inhibit corrosion, or whe- 
ther passivity will practically not occur because of a ra- 
pid dissolution rate of the surface film. This inconve- 
nienco is, however, eliminated by means of the polariza- 
tion curves: electrode potentials vs. current densities 
applied to the electrode. Therefore it is appropriate to 
consider that the static aspect stresses the potential-pO*— 
diagrams and that the dynamic aspect stresses the po- 
tential-current density curves. 

The construction of the equilibrium diagrams E-pO*— 
is based on the well-known method used by Pourbaix in 
his E-pH diagrams (Pourbaix 1949). The principle of the 
graphical representation E-pO*— is quite simple, but its 
application sustains a certain number of difficulties. In 
fact, if one wants to draw the diagram for a metal-oxya- 
nion melt system, for example, it is necessary to outline 
the complete inventory of all the likely products of for- 
mation and to know the thermodynamic data relevant to 
these compounds. This basic work being done, it will be 
necessary to examine all the possible reactions between 
the compounds considered. In order to make the cons- 
truction as systematic as possible, any reaction in which 
substance À is changed into substance B is written in a 
form such that the only other species appearing are sol- 
vent ions, R*, oxide ions, O*—, and anions, e—. Such 
reactions therefore have the following general form, all 
symbols being collected on one side: 


aA+bB+cR-mO'- +ne-=0 (26) 


end it was shown in section on The single potential of an 
electrode that the corresponding conditions of equilibrium 
can be expressed in the form 


* The symbol R is used here to represent the reducilbe or and oxidizabla species resulting from the acid-base dissociaton of the 


ruelt. 
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aua + bug + C4pT MEG — n F E (06 =0 


E 0,96 being defined by the relation (7), which for the 
present case becomes: 


23 RT 
oo = East ——p [2 log tA) +b log (8), 
2.3 RT 
+ € log 1R)] + - E m pO*— (27) 


It is therefore possible to consider several classes of 
reactions: 


1. Reactions in which e” anions do not take part. 
Equilibrium in these reactions depends or does 
not depend on pO*— depending on whether m 
+ O or m = O respectively. These reactions are, 
of course, purely chemical reactions. 


2. Reactions in which e” anions take part. 

These reactions are electronation-deelectronation 
reactions; therefore they are electrochemical 
reactions. Equilibrium in these reactions is given 
by expression (27). If the coefficient m is equal 
to zero the equilibrium condition (27) is inde- 
pendent of pO'—; if m-+ o equilibrium depends 
on both pOº— and potential. 


in these conditions, on the diagram one meets a series 
of staight lines having several. positions. For example, 
reactions without de-electronation involving O*— ions are 
represented by vertical lines; reactions with de-electrona- 
tion in which O*— ions do not take part are represented 
by horizontal lines; in the other cases, the lines will have 
variable slopes depending on the values m and n. These 
group of lines divide the field of the diagram into several 
domains, each corresponding to the domain of thermo- 
dynamic stability of the particular species. 

Let us reproduce as an example the essential steps 
relative to the construction of a diagram E/pO*— for Ni 
in Na.SO, at 900ºC. Following the above principles, the 


TABLE | 


Frea energy data for the species present in the Ni — Na; 
SO, system at 1173K 


Formula Standard State —uº (kcal/mole) 
Ni solid metal 13.5538 
N+ Solid metal ion 56.4713 
Ni O solid oxide 75 8961 
Nis solid sulphide 47.4616 
NiiS; liquid sulphide 118.0500 
O, gas, 1 atm. 63.033 
Ore solid oxide ion — 74.1327 


e solid sulphide ion — 88.3458 
SO. gas, 1 atm. 148.4475 
SO, gas, 1 atm. 177.5087 


ss 


energetically probable products of the assumed reactions 
between the melt and the metallic nickel in contact with 
it aro envisaged. The species and their free energies of 
formation are given in Table |. The standard free energies 
of formation (u? )at 900ºC for the species considered 
were calculated from equations given by Nagamori and 
Ingraham (1970) and Kubaschewski et al. (1967). The 
question of existence or non existence of the compounds 
considered in the table is not relevant here, because the 
example is only presented to illustrate the procedure for 
constructing the diagrams. Details of this particular case 
are given by Sequeira and Hocking (1975). Therefore 
the reactions are classified following the adopted princi- 
ple and were assembled in Table Il. 

Figure 4 shows, the conditions of stability for nickel, 
its oxides, sulphides and ions, and the domain of ther- 
modynamic stability of the sulphate melt at 1 atm. total 
pressure (assuming that the three extreme lines (b), (c), 
(d), for which P5, = 1atm., Peo,= tatm. and log 


(S—)| = —6 respectively, determine the potentials limi- 
ting the stability of the sulphate melt). To make easy the 
understanding of these diagrams the lines referring to the 
several reactions will carry the same letter as the equi- 
librium equation they represent. Moreover, when neces- 
sary, the families of lines referring to a particular equi- 
librium where the ionic activities have certain given va- 
lues will be drawn. For example, for the reaction (a), one 
meets the family of lines referring to log (Ni'+] values of 
O, —2, —4, —6. In this case, activities are indicated on 
each line of the family. If it is assumed, as by Pourbaix 
(1949) for aqueous solutions that a metal does not cor- 
rode unless the substance forming on its surface has a 
solubility in excess of an arbitrary * value of 10-* in the 


TABLE 1 


Electrochemical reactions and their equilibrium con- 
ditions for the Ni — Na.SO, system at 1173 K 


Reactions in which e anions do not take part: 
(a) NiO = Njt + 0O'— 
log (N+) = — 17,33 + pO— 


Reactions in which e— anions take part: 


(b) 50,+26-=0- 
E, = — 2.2891 + 0.0582 log P 9 + 0.1163 pOº 
(c) SO. +- 26 = SO. + e 


E, = — 4.5135 — 0,1163 log Psoo + 0.2326 po*— 


* This limiting concentration is suggested by the argument that from initial corroslon a concentration of 10 “ of metal lons 5 es- 
tablished very rapidly in the electrolyte. If the electrolyte is not connuously renewed, the maintenance of this metal concentration does 


not represent a serious loss. 
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Fig.5 — Theoretca|l conditions of corrosion, immunitv and passivat ion for nickel 
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(d) SO, + Be- = S'- + 30 

E, = — 3.2145 — 0.0291 log (S—) + 0.1163 pO- 
(e) NIO +2e- =Ni+0O>- 

E, = — 2.9574 + 0.1163 po 
(f) 3Ni+2S0,+126- =NiS: +60 

E, = — 3.3677 + 0.1551 po 
(g) NiS: + SO, + 6e-=3NIS +30 

E =— 34714 + 0.1551 po- 
(h) NE+:SO,+8e- =NiS+30- 

E, = — 2.7848 + 0.0291 log (Ni*t] + 0.1163 po 
(i) 3Ni+ + 2S0, + 18e- = NiS. + 6 0 

= — 2.5551 + 0.03876 log (Ni*+) + 0.10337 pO*- 
(j) 2NiO + 250, + 186 — = NiS. + 90% 

E, = — 3.2309 + 0.14214 po”: 


appropriate concentration units, the line representing this 
solubility on the potential — pO* diagram divides it into 
two dom” ins; one of non-corrosion and one of corrosion. 
Ye: similarly to aqueous systems (Pourbaix 1949), the 
domains of noncorrosion can be divided into domains 
oi immunity (where the surface of the metal only con- 
sists of the metal itself) and domains of passivity (where 
the metal is covered with a complex film). The existence 
of tho domains of stability shown in Fig. 4 and the 
assumptions above allows the prediction that, accor- 
ding to the po” of the melt and the electrodo 
potential of the nickel, it will be, as shown in Fig. 5 
(which is easily derived from Fig. 4), in a state 
of total protection (by either immunity or perfect passi- 
vation) or of localized corrosion (in the case of imperfect 
passivation) or of generalized corrosion (in the domain of 
corrosion). 


Further details on these diagrams are beyond the 
scope of this paper and for the practical employment of 
these methods the reader is referred to the several works 
published over the last few years (Baudo and Tamba 
1969, Bombara et al. 1968, Busson et al. 1967, Conte 
and Ingram 1968, Lewis 1971, Littlewood 1962 a, b, 
Rahmel 1968, Sequeira and Hocking 1975, et al.). 


It should be noted once more that the domains of 
thermodynamic stability considered above give only a 
theoretical possibility of oxistence, not a certainty. In par- 
ticular as far as passivation is concerned, it must be 
pointed out that a proper passivity in molten media has 
not been obtained up to date, perhaps due to special 
properties of the film involved, mainly adherence and 
coherence. 


POLARIZATION MEASUREMENTS AND FREE CORRO- 
SION RATES 


In this section an account is given of the available 


methods for obtaining the potential-current curves as 
well as the free corrosion rates. 
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It is customary to classify the electrochemical me- 
thods for obtaining such potential-current curves into 
potential-sweep, current-sweep, potential-step and cur- 
rent-step methods. The former method, sometimes called 
potentiodynamic (or potentiokinetic) polarization, is ac- 
complished by continuously changing the applied elec- 
trode potential (with feedback control ) at a constant rate 
and simultaneously recording current. In the current- 
-sweep (intensiokinetic) methode, the electrolysis current 
(with feedback control) is varied with time according to 
a given programme, and the electrode potential is re- 
corded as a function of current. In the potential-step 
method, electrode potential is rapidly changed over a 
finito increment, and current measured after a predeter- 
mined time interval, and this process is repeated. In the 
last (above) method, a current step is applied, and the 
resulting variation of potential with time is recorded. 
These two step methods are also referred to as potentios- 
tatic and galvanostatic (or intensiostatic) respectively, to 
indicate that at any instant E (or |) is held constant, by 
means of feedback control. 


All the polarization methods described above are 
well established in the ambient temperature range (Ede- 
leanu 1958, Greene 1962, Greene and Leonard 1964, 
Pourbaix and Vandervelden 1965, and Prazak 1963) for 
electrochemical studies of metals, the best or most ac- 
curate of them depending of course, on its intended use: 
if the purpose is to study kinetics phenomena, the me- 
thode in which the electrode potential is controlled are 
the best; when the permanent physico-chemical modifica- 
tions, due to the applied current, are the objective, the 
methods in which the electrolysis-current is controlled are 
preferable. Generally speaking, however, potential-step 
and potential sweep methods will be more useful for 
study of corrosion phenomena than current-step and 
curren-sweep methods, because they give more precise 
indications concerning the quality of the passivation and 
concerning the conditions of activation. The investigation 
of corrosion processes in melts by these standard electro- 
chemical techniques seems more limited (Arvia et al. 
1971, 1972 a, b, Baudo et al. 1970, Casino et al. 1971, 
Davis and Kinnibrugh 1970, Delarue 1959, Janz and 
Conta 1964, Kazantsw et al. 1968, Pizzini and Agace 
1965, Pizzini et al. 1964) than in aqueous systems; how- 
wever, almost all workers in the field of high-temperature 
agree that the polarization methods appear to be quite 
suitable in corrosion studies at high temperatures in 
molten electrolytes. 

The polarization curves have several shapes, but 
near the potential! axis their character is similar to that 
shown by the curves plotted in Figures 1 and 2. 


A potentiodynamic polarization curve provides a gra- 
phic summary of the corrosion characteristics of tested 
material so that its plot and evaluation can be used as a 
very rapid and practically nondestructive corrosion test. 
However, the corrosion characteristics obtained are ac- 
curate only for the medium in which the curve was. 
obtained. The importance of the method, therefore, does 
not lie in any final evaluation of the rates of corrosion in 
a particular medium, but rather in the rapid evaluation of 
fundamental corrosion properties with the possibility of 


sensitive, relativa comparison (Prazak 1963). Another 
featura of the method is that analysis of such potential 
-current curves, together with the E-pO— diagrams, ena- 
bles us to predetermine those reactions theoretically pos- 
sible and impossible for each electrode potential value. 


One of the most important electrochemical values 
for the determination of the kinetics and mechanisms of 
primary corrosion processes is the value of the free 
corrosion current. When the pimary corrosion process is 
slow (absence of concentration polarization) the free cor- 
rosion current can be obtained by potentiodynamically 
polarizing the alloy electrode at a very low sweep rate, 
in the neighbourhood of the free corrosion potential. The 
alloy-dissolution rate is obtained by plotting potential vs 
log current and extrapolating the linear anodic and catho- 
dic branches to the free corrosion potential, the intercept 
giving the free corrosion current. This method can also be 
applied when some (relatively small) concentration pola- 
rization is present; the latter can be taken into account 
by using the approximate expression: 


| * yo 


meas, 
| 


| d  meas, 


Es 
— 


whero legs is the measured current in the presence 
of concentration polarization and |, is the diffusion cur- 
rent (of course this correction implies the exact know- 
ledge of the diffusion current). When the primary corro- 
sion processes are sufficiently fast, relaxation methods 
must be used. The objective of these methods is the 
study of the electrode processes before they become dif- 
fusion controlled, Tha experimental techniques associated 
with these methods are very complex, but they are the 
most important ways which permit accurate determina- 
tions of fast charge transfer corrosion rates. The princi- 
pal relaxation methods for these measurements have been 
considerably revised in the last few years (Graves and 
Inman 1970), and a useful exposition of the «classical» 

relaxation methods is that of Damaskin (1967). The 
selection of the method to measure the free corrosion 
currents, depending precisely on rate phenomena, is ob- 
viously difficult (because the values of these corrosion 
rates are not known a priori). 


In general the high temperatures of fused salt sys- 
tems cause most corrosion processes to be rapid (rever- 
sibla with respect to the non-relaxation methods). Gene- 
rally speaking, thereby, relaxation methods will be the 
most appropriate for the determination of corrosion rates 
in melts. 


The measurement of free corrosion rates in molten 
salt systems is reported by Baudo et al, (1970), and by 
Cutler (1971). The first workers determined the corro- 
sion rates using the conventional method, extensively and 
successfully appkied at low temperatures in aqueous sys- 
tems; they presented the results in the form of Tafel 
plots but they do not give information on the current- 
-potential measured values from which the plots of E vs 
log i were obtained. It is rather surprising in view of the 
relatively high corrosion rates quoted that it had been 
possible to obtain pure Tafel zones. It is certain that this 


method enables the measurement of rough values (more 
of qualitative than quantitative value) of the free corro- 
sion rates, but it is adviseable to put in practice LaQue s 
opinion that new and more accurate tests should be intro- 
duced for a better understanding of the corrosion mecha- 
nisms (LaQue 1969). A cyclic potential sweep method 
is used by Cutler (1971) to determine the free corrosion 
rates which seems scientifically more appropriate than 
the conventional one. 


The search for an accurate method for the determi- 
nation of the kinetics and mechanisms of corrosion pro- 
cesses in molten salt systems has led to a study of the 
recent «aqueous» literature, from where the double pulse 
potentiostatic (D.P.P.) method improved by Barbier et al. 
(Amosse et al, 1968, 1970, and Bouteillon et al. 1968) 
was selected. This method enables a direct reading of 
the charge-transfer current, it being possible to establish 
a linear relationship between the potential and the loga- 
rithm of the current (which is especially good for extra- 
polation of the free corrosion currents). Theoretical de- 
rivations as well as experimental verifications of the 
merits of the D.P.P. technique to the measurement of the 
free corrosion rates in melts will be provided by a for- 
thcoming paper. 


CONCLUSION 


In this paper, a glance at the essential elements con- 
cerning the electrochemistry of corrosion phenomena in 
molten salts was given, The electrochemical equilibrium 

diagrams as well as the potential-current curves were 
shown to be two important tools put at the disposal of 
the scientist for the solution of corrosion problems. These 
diagrams, which are of a thermodynamic character, can 
serve primarily as a guide in carrying out such polariza- 
tion curves and in the interpretation of the results ob- 
tained. 

The techniques which can be used to obtain the po- 
larization curves as well as the corrosion rates were brie- 
fly described. lt was shown that the potentiodynamic 
measurements are particularly important in the evaluation 
of a series of experimental alloys. As far as passivation 
is concerned, analysis of such measurements, together 
whith non-electrochemical examinations of the corres- 
ponding surfaces (e.g. by means of visual observa- 
tion, electron probe micro-analysis) should inform on 
their true passivation. 


The study of the electrode parameters (exchange 
currents, Tafel slopes, and equilibrium potentials of the 
metal-dissolution and electronation reactions) in melts 
requires the use of the recent relaxation methods; in 
view of the dependence of the free corrosion rate on 
these parameters (as is clearly brought out in Figure 3), 
its accurate measurement must be made using such te- 
chniques. 
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SYMBOLS 

A Area of electrode 

a Constant 

b Tafel coefficient, Tafel slope (0 E/ d log) 

E Reactional electrode potential or electric potential 
of an electrode when it carries a current den- 
sity (i) 

É. Free corrosion potential 


Em Mixed potential 

Single potential of an electrode or electrical equi- 
librium potential, measured with respect to a 
standard reference electrode 

EQ Normal electrode potential, or equilibrium potential 


E, when the reactants are all in the standard 


state 

em Electrons 

F Faraday's constant (96.5 KC/mol or 23.06 Kcal/ 
/volt-mole) 


A G * Chemical free energy of activation 

a G* Electrochemical free energy of activation 

| Total current (i x area) 

| Free corrosion current 

Diffusion current 

i Any reacting species; current density 

Exchange current density 

Equilibrium constant 

M;, Chemical symbol of reactant i 

1 M | Activity of constituent i (thermodynamic concen- 


tration) 

n Valency of metal, number of electronic charges on 
solved-cation 

P Pressure (in atmospheres) 

pO*- Cologarithm of oxygen-ion activity 

R Resistance of the electrolyte; gas constant (1.98 
cal/“C) 

T Absolute temperature 

a Symmetry factor 

n Overpotential 

H; Chemical potential of constituent i 

u Electrochemical potential of constituent i 

nº Standard chemical potential of constituent i 

p Inner potential 

rp; Stoichiometric coefficient of constituent i 


Subscripts such as in Ein indicate that the 
symbol refers to reaction n 
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O modo como o homem adquire conhecimento do universo, 
como o conceptualisa e actua. 
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UTILIZAÇÃO DE MODELOS AERODINAMICOS NO ESTUDO 
DE PROBLEMAS FLUVIAIS 


Técnica N.º 429 — L — 10-1975, p. 13 a 24 


A utilização de modelos asrodinâmicos para o estudo de 
escoamentos com superfície livre é uma técnica que parece 
ser digna de interesse, por permitir estudar certos proble- 
mas fluviais em modelos de muito menores dimensões do 
quo us dos modelos hidráulicos tradicionais proporcionando, 
assim, grandes economias. Por isso, foi decidido empreender 
no LNEC, o estudo dos aspectos teóricos e práticos do es- 
tabelecimento e exploração dos modelos aerodinâmicos, 
com vista à avaliação das suas potencialidades e, eventual- 
mente, à preparação da sua futura utilização de rotina. 
Neste artigo faz-se um breve resumo histórico dos modelos 
aerodinâmicos, refere-se o equipamento, materiais e técnicas 
de construção a que se recorre, apresentam-se os critérios 
de semelhança mais significativos, dá-se ideia dos con- 
dicionamentos e aplicabilidade desses modelos, descrevem- 
-So as operações de taragem e exploração e comparam-se 
as vantagens e inconvenientes relativos dos modelos asro- 
dinâmicos e dos modelos hidráulicos tradicionais. Por úl- 
timo relata-se a actividade desenvolvida e em curso no 
LNEC no âmbito dos modelos aerodinâmicos. 
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SOBRE A APLICABILIDADE DUM MODELO 
A PREVISÃO DA POLUIÇÃO ATMOSFÉRICA 


Técnica N.º 429 — L -—10-1975, p, 25 à 32 


Esto artigo destina-se a pessoas que, não estando muito 
familiarizadas com modelos de previsão de poluição atmos- 
férica, estejam mais interessadas em conhecer as suas po- 
tencialidades do que as bases físicas do processo. 
Descreve em termos gerais as utilizações mais correntes 
destes modelos nomeadamente planificação da utilização do 
solo, avaliação de estratégias de controle e funcionamento 
de sistemas de alerta. 

Os modelos Gaussianos, apesar das suas limitações, têm 
permitido obter resultados satisfatórios, são de utilização 
fácil o necessitam pouco tempo de computador, pelo que 
lhes é dado especial relevo. Referem-se os dados necessá- 
rios &o funcionamento destes modelos assim como as in- 
formações que podem fornecer. 

No N.E.E.M. foi desenvolvido um modelo Gaussiano de 
fontes múltiplas e tem sido realizado trabalho conducente a 
um modelo geral de previsão de poluição baseado na reso- 
lução das equações da camada limite atmosférica e de di- 
fusão turbulenta. 
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BASES ELECTROQUIMICAS DA CORROSÃO METÁLICA 
EM MISTURAS FUNDIDAS 


Técnica N.º 429 — L — 10-1975, p. 33 a 46 


Uma vez apontada a natureza electroquímica da corrosão 
metálica em meio fundido, descrevem-se os parâmetros 
teóricos essenciais da termodinâmica e cinética electroqui- 
mica da corrosão. Dá-se relevância à electrocinética da 
corrosão; com base na teoria de «Erdey-Gruz-Volmer» e nos 
conceitos de ácido e base segundo Lux-Flood, estabelece-se 
uma relação entre o potencial de corrosão e o cologaritmo 
da actividado exigeniónica do electrólito fundido. Introdu- 


zem-se os diagramas de equilíbrio E-pO?- , cuja constru- 
ção se faz segundo o método de Pourbaix, de representação 
do dados termodinâmicos; exemplifica-se com a construção 


de um diagrama E-pO'” para o sistema Ni-Na,SO, a 
900º C. Finalmente, sugere-se que se encare a corrosão me- 
tálica nas misturas fundidas, em termos da cinética electro- 
química, recorrendo-se às técnicas transientes, usadas em 
electroanálise, para o estudo da cinética e mecanismos dos 
processos corrosivos. 


Synopsis of articles published in «Técnica» n.º 429 
L — Qrnonar id 


U. D. C. 628.5 
M. TERESA PONTES 


A NOTE ON THE APLICABILITY OF A MATHEMATICAL 
MODEL TO THE PREDICTION OF ATMOSPHERIC 
POLLUTION 


Técnica Nº 429 — L — 10-1975, p. 25 a 32 


This paper is mostly intended to potential users of air pol- 
lution models and to those more interested in solving a par- 
ticular problem than on the physical processes involved. 
It describes in broad terms the more common use of pol- 
lution models namely planification of land use, survey on 
control strategies and alert systems. 

Gaussian models are particularly valuable for solving this 
class of problems as, in spite of their limitations, they lead 
to adequate results, are simple to use and require little 
computer time. Reference is made to the input data requi- 
rements and to the type of information the model can pro- 
vide. 

A multiple source Gaussian model was developed at NEEM. 
Work on a general model based on the solution of the 
atmospheric boundary layer and turbulent diffusion equaiions 
is under way. 
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«ELECTROCHEMICAL BASIS OF METALLIC 
CORROSION IN MELTS» 


Técnica N.º 429 — L — 10-1975, p. 33 a 46 


After pointing out the electrochemical nature of corrosion 
reactions in molten salts, am account of the essential para- 
meters concerning the electrochemical corrosion phenomena 
is presented. Emphasis is given to the electrodics of cor- 
rosion; based on the «Erdey-Gruz-Volmer» theory and on the 
«acid-base concept» of Lux and Flood a relationship bet- 


ween the corrosion potential and the function pO “is deri- 
ved. Diagrams in which the equilibrium potential E in the 
metal-fused salt system is plotted against the activity of 
oxide ion in the melt are developed, using the Pourbaix 


method of data presentation; as an example, an E-p0'” 
diagram is constructed for nickel at 900º C, in contact with 
molten sodium sulphate, and its limited nature as an aid 
to corrosion control is summarized, The need for study of 
metal corrosion in terms of electrochemical kinetics is in- 
dicated, Further, it is suggested to apply the newer electro- 
chemical tecchniques to the study of corrosion kinetics in 
high-temperature media. 
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How man aquires knowladge about the universe concopiua- 
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The uso of air-flow models for studving open-channel flows 
seems an interesting techniques, because it permits the 
study of some river problems with smaller models than the 
traditional hydraulic ones, thus affording great economy. 
So, it was decided to start at the «Laboratório Nacional de 
Engenharia Civil» the studv of theoretical and pratical as- 
pects of establishina and using air-flow models, to evaluate 
their potentialities and, eventually, to prepara their future 
routine use. 

The paper makes a brief historical summary of air-flow 
models, refers to equipment, material and construcion 
techniques, shows the most important similar criteria, gives 
an idea of the restrictions and applicability of air-flow mo- 
dels, describes the problems of its calibration an use and 
compares tha advantages and disadvantanes of air-flow mo- 
dels with thosa of traditional hydraulic ones. Finally it re- 
ports on the activity which has been carried out and that 
is in progress at the LNEC in that field, 
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